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1. DOTYCHCZASOWY STAN BADAŃ NAD EFEKTAMI TOPOGRAFII 
POWIERZCHNI HETEROGENICZNYCH 


Mija już prawie pół wieku od momentu, gdy Langmuir 
zaproponował swoje słynne równanie izotermy adsorpe ji/ 7, 
które do dziś nie straciło użyteczności w badaniach 
zjawisk adsorpcyjnych. Mało kto dziś jednak pamięta, 
ze był on też pierwszym, który zwrócił uwage na ważne 
efekty, jakie moze powodować niejednorodność energetyczna 
powierzchni. Swoje równanie uogólnił na przypadek , gdy 
na powierzchni adsorbentu można wyróżnić kilka rodzajów 
centrów adsorpcyjnych, charakteryzujących się różnymi 
energiami adsorpcji. 

Ta część pracy Langmuira poszła w zapomnienie w latach 
trzydziestych, bowiem ówczesny trend naukowy preferował 
sprawy oddziaływań między molekułami zaadsorbowanymi e 
Wówczas to sformułowano bardzo dziś popularne równanie 
Fowlera-Guggenheima dla monowarstwowej adsorpcji 
zlokalizowanej z oddziaływaniami typu "najbliżsi sąsiedzi", 
Już wkrótce, bo na początku lat czterdziestych sformułowano 
nie mniej popularne równanie Hilla-de Boera/ */ dla 
monowarstwowej adsorpcji mobilnej z oddziaływaniami 

w fazie adsorpcyjnej. Jednakże mimo dalszego doskonalenia 
Opisu teoretycznego efektów oddziaływań między molekułami 
zaadsorbowanymi, notowano coraz to liczniejsze fakty 
doświadczalne, których modele adsorpcji z oddziaływaniami 
nie były w stanie wytłumaczyć. 


Dla przykładu, powszechnie obserwowano niezrozumiały, 
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bo szybki wzrost adsorpcji w obszarze niskich ciónier ` 

i bardzo powolny dalszy wzrost adsorpcji w obszarze 
średnich i wyższych ciśnień. Bardzo rzadko eksperymentalne 
izotermy miały skoki, które możnaby interpretować jako: 
przejścia fazowe przewidywane przez teorię, np. gaz 
dwuwymiarowy-ciecz dwuwymiarowa, jak w równaniu 

Hilla-de Boera. 

Kalorymetryczne pomiary izosterycznych ciepeł adsorpcji 
dawały zazwyczaj ciepła malejące wraz ze wzrostem pokrycia, 
czyli wynik wręcz odwrotny do przewidywań teorii bazującej 
na modelach adsorpcji z oddziaływaniami. Pomiary pojemności 
cieplnych filmów adsorpcyjnych nie wykazywały nieciągłości 
przewidywanych przez t4 teorię w obszarze niskich 
temperatur i wysokich pokryć powierzchni; dla niższych 
pokryć powierzchni obserwowano na ogół duże odchylenia 

od zachowania przewidywanego przez przybliżenie Einsteina 
czy Debye’a. Prace eksperymentalne na ten temat można 
byłoby cytować dziesiątkami, lecz to nie wchodzi w zakres 
tej pracy. 

Gdyby powyższe fakty określić jako zupełnie niezrozumiałe 
W owym czasie, byłoby to stwierdzenie nieprecyzy jne. 

Bo oto, niektóre z nich można było zinterpretować 

w oparciu o koncepcję powierzchni heterogenicznych, uje ty, 
matematycznie w pracach E NSF 1 ETN roza” W 

z lat trzydziestych. Przy dość słabej wymianie informac ji 
naukowej w tych latach, prace tych radzieckich uczonych 
nie odbiły się należytym echem w ówczesnym ówiecie naukowym. 


Na gruncie radzieckim zapoczątkowały one świetny okres 
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rozwoju badań w dziedzinie adsorpcji, których wyjątkowa 
owocność przypada na lata ezterdzieste. Zasługą wspomnianych 


badaczy było między innymi przyjęcie matematyczne j koncepcji 
rozkładu energii adsorpcji w -— funkcji ciugtej. 
Pozwolizo to dokonaé prawdziwego skoku w formaligmie opisu 
zjawiska adsorpcji na powierzchniach heterogenicznych. 
To właśnie Todes 2 przekształcił po raz pierwszy röwnanie 
całkowe z lokalną izoterma Langmuira do postaci 
transformaty Stieltjesa. Z powodu nieznajomości tej pracy 
opublikowane j w czasopiśmie ukraińskim, wyniki jego 
w trzy lata później powtórzył Sipa/ ۰ 

W latach czterdziestych i na początku lat pięśdzie- 
siątych mnożą się prace oparte na całkowym równaniu 
dla globalnej izotermy adsorpcji, lecz mimo to waga 
problemu nie zawsze jest należycie rozumiana, Bo oto sam 
Todes// a również Honig i Rosenbloon/ D uważają; | 
że globalna izoterma adsorpcji jest zupełnie nieczuła 
na postać funkcji rozkładu energii adsorpcji. 
W latach sześdziesiątych rola heterogeniczności powierzchni 
w zjawiskach adsorpcyjnych nie budzi żadnych wątpliwości. 
Wybitni uczeni w tej dziedzinie zwracają uwagę na fakt, 
jak niewiele jeszcze zrobiono na tym polu. | 
W dyskusji z profesorem Kisielowem, Everett! Y uvaza, 
ze "przed nami jest jeszcze długa droga do zrozumienia 
adsorpcji na powierzchniach heterogenicznych". Byto to 
اس و ین‎ lat temu, a w sześć lat potem Pierotti w swoim 
artykule ل‎ ol a stwierdza, ze "problem hetero- 


genicznoóci vozostaje wciąż jednym 2 wielkich nierozwiązanych 


wi Lifted wén vz -sqiecke (çan ağa $ 


"> TORO a mil 


ne.” 


problemów adsorpcji". Mimo niewatpliwych postqpów 
na tym polu, jeszcze cztery lata temu Sams? ۷ uważał, 
że "problem stworzenia zadawalających teoretycznych 
podstaw dla opisu adsorpcji na powierzchniach hetero- 
genicznych jest niewytpliwie jednym z najbardziej trudnych 
1 skomplikowanych problemów w adsorpcji fizycznej". 
Złożoność problemu i mały stopień jego znajomości 
Sprzy jaz też powstawaniu poglądów ekstremistycznych 
w pewnym sensie. Ulegali im nawet badacze tej miary jak 
Adamson! 1*/ , który uważał, ze "dla większości powierzchni 
adsorbentów rozkład energii adsorpcji jest charakterystyką 
dominującą, jednoznaczną, i na nią należy kłaść główny 
nacisk w opisie teoretycznym adsorpcji, niż np: na model 
adsorpcji lokalnej". RE ə wyznawał zupełnie skrajny 
poględ: "postać przyjętej izotermy lokalnej nie ma 
żadnego wpływu na postać wyznaczone j funkcji rozkładu 
energii adsorpcji". Godny zanotowania wstrzemięź liwość 


Stek Mé 7 WH 
1 intuicję reprezentował Halsey 


który jeszcze‏ و 

w 1952 roku pisał: "Adsorpcja fizyczna reprezentuje rodzaj 
mariage de convenance, ponieważ zarówno oddziaływania 

w fazie zaadsorbowanej jak i heterogeniczność powierzchni 
są niezbędne do zinterpretowania typowych izoterm", 
Wzięcie obydwu tych czynników pod uwagę, stawia nowy, 
niezwykle PA problem wpływu topografii powierzchni 

w adsorpcji na powierzchniach heterogenicznych. Chodzi tu 
Age o to czy istniejy jakieś korelacje (w sensie 
rozmieszczenia na powierzchni miejsce adsorpcy jnych) 


między centrami adsorpcy jaymi o zbliżonej energii adsorpcji» 
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Mozna tu wyróżnić dwa skrajne przypadki: 
W pierwszym przypadku, centra o jednakowej energii 
zgrupowane są w "płatach"; a więc istnieje pełna korelacja 
w rozmieszczeniu centrów adsorpcyjnych. Płaty te są 
na tyle duże aby w opisie każdego z nich móc stosować 
metody termodynamiki statystycznej oraz by można było 
zaniedbać oddziaływania pomiędzy molekułami zaadsorbowanami' 
na różnych płatach. W obrazie takim, faza adsorpcy jna | 
(układ molekuł) składa się z dużej ilości niezależnych | 
od siebie podukładów. Uważa się jednak dalej, że w sensie 
praktycznym podukłady te sa na tyle jeszcze maze, by widmo 
dyskretne można było nadal zastępować rozkładem energii | 
adsorpcji w postaci funkcji ciągłej. Model płatowy został 
po:raz pierwszy sformułowany jasno przez T Sarina i ۱ 
Lewicza/ 1>/ , oraz niezależnie, przez Halsey"a i Taylora/ 19 . 
Najbardziej wyczerpujace studia nad tym mode lem topogra- ` 1 
ficznym wydaja sie 6 dzierem Rossa i Oh وی‎ 17/7 
Oni tez byli autorami przymiotnika "homottatic" dla tego 
modelu, choó w literaturze przyjęło się raczej określenie 
"patchwise". | 
Drugim, krańcowo różnym modelem powierzchni hetero- 
genicznej jest założenie o zupełnym braku korelacji 
w rozmieszczeniu centrów adsorpcyjaych o podobnej energii 
adsorpcji. Innymi słowy, uważa się, ze centra o różne j 
energii adsorpcji rozmieszczone są na powierzchni 
w sposób zupełnie przypadkowy. W takim obrazie układ 
molekuł zaadsorbowanych musi być traktowany jak jedna =< 


termodynamiczna całość. Taki model określa sie w litera- 
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-turze jako topografie typu "random", 

Jest w istocie rzeczą zastanawiajaca, dlaczego ponad 
90% prac dotyczących adsorpcji na powierzchniach ne take. 
genicznych zakiada istnienie topografii typu piatowego. | 
W gruncie rzeczy nie ma zadnych dowodów eksperymentalnych 
na to, ze ten typ topografii jest ciii ce SE dla 
większości układów adsorpcyjnych z powierzchnią hetero- 
geniczną. Pewnym usprawiedliwieniem mogą tu być racje 
historyczne; gdyż model ten został sformułowany jako 
pierwszy, co dawało mu pewien "handicap" w popularności 
wśród badaczy. Natomiast, model powierzchni o rozkładzie 
przypadkowym miejsc adsorpcyjnych wchodził bardzo powoli 
do literatury. 
Rzecz raczej dziwna, bo jeszcze w latach czterdziestych 
Hi11/ 18/ uogólnił przybliżenie Bragza-Williamsa 
i Bette-Guzgenheima (quasichemiczne), na przypadek 
powierzchni heterogenicznych z topografią typu "random". 
Na początku lat pieddziesiatych fómpkins/ "W rozważał 
zachowanie się dwuwymiarowego gazu Van der Waalsa 
na powierzchniach typu "random", lecz były to poczynania 
zgoła odosobnione. 
Dopiero Steele/2V w 1965 roku podniósł zdecydowanie 
problem topografii powierzchni stwierdzając: "ścisłe 
określenie własności układów z heterozeniczną powierzchnią 
wymaga kompletnej znajomości topografii całej powierzchni". 


Steele był też pierwszym, który z danych eksperymentalnych 


dotyczących trzeciego współczynnika wiriału ¿az-cialo Stare, 


próbował wyciągnąć określone wnioski na temat topografii 


powierzchni. Jego zdaniem, założenie topografii typu 


"۲ N 


"A 
Ta 
Ly [ , Y; " 1 di 
“x = Iz Ë ہے‎ de 
0 ن‎ M HS Dy | 
< ۳ھ‎ 9 Ag 
«we P rØ £ Lait AP 3 یی | سد‎ 
= Í 5 1 J vty | k. M 
EX = P 1 5 » T ` - à 
B 5 . 4 " f Psy ge? ; ۸ 
j JA is (7 Lü % ~ d éi 
` 
- - +4 
174 1 ) GEI) [M i J Bs i 1 ابا‎ 7 
l | Wot xoqaoaba wobafXs Lohoraxplw 
` kd ۱ )4 t Il 一 AN D 2 e Eu isn 
A DAN 1 Ly Ton rz pr n wf. r ura KL A td 
A 
* pə اس‎ + < l - 4 bn | 
؟‎ 1 7 g Dot zg ul nun OF وس 1 ہہ‎ ٤ 
= =Ë E a š k A db 4 4 a ' : m. > e ۱ 
1 ۳ د‎ 1 9۵0۵۵ Ipa LATTAN -¥soaher DOIBW ! 
a d w .4 ion a ١ Pr a Y 1 w ۳ «ge 
ä ۹ te ارزع هبني‎ Aaaa mca انا‎ ta 


UST" əyan ce of 


U 


00 : ly „r T | E yv 9b 2 "ame ts L 8 2ھ‎ 
SITE بس‎ Sil fies من ری‎ “S” WISH 

beg لکن‎ og ¿Lancia data) sari uranın P3 
aan uryf ge s تاهو دمح + اج م‎ arl ¿odos eitwog ? 
ŻE Lo m ty : ECT» DE "mni: وو جم‎ n raj müs ap? 

- W) 


Y ün qul ken ağanı re gla 5 


pinag uw kond, o Boe qon ere woz Sp | 


تید rəb bor‏ ج1 cost‏ عبج - مل اتھ۔ 
a ١ | = vu Q E‏ سک - 
ARE i |‏ سے LINE‏ ۱ 
vip kət = Vis, f ansa. osa Saq |‏ ریہ نا E "fat aşya bad‏ 


ved pta" LM ln dn & A پوس ہتفہ‎ ku yoqa 20 


: ۴ pin Dah) ELDER SENT ` RODE “Azra İv eiculba | 
b -—— و9‎ reus KZ 3 uż dot «tle L*^o25igmox “yu ۱ 
. . Oy feodi Coe ai a xan 2.5... die 
2—- odu in. بیش‎ SCHT ۱ uin و مو شش سے‎ ACA | 
1 lg is Pt «© ۸۰۵۵ calado ao Soo dig p 
- AT. ۶ —. Sitti Slam ağladar dari 

7 "gra Y JM oc = <4J BEM 5 


1 ui 


alqı 


"random" prowadzizo, przynajmniej w przypadku niektörych 
8dsorbentów węglowych, do lepszej zgodnoóci teoretycznych 
i eksperymentalnych wartoáci 3- -go współczynnika wiriału 
gaz-ciało stałe, | = 

Badania tego typu kontynuowa i następnie Pierotti 
1 Thomas” 217 , Rudziński/ £7 SS Rippa i Zgrablicn/ 29 
byli pierwszymi, którzy na gruncie formalizmu wirialnego 
rozpoczęli rozważania nad pośrednim modelem topograficznym, 

modelem, w którym zakłada się istnienie częściowych 

| jedynie korelacji położeniowych między centrami 
o jednakowej energii adsorpcji. 
Jak z tego zestawienia wynika, liczba prac dotyczących 
efektów topografii powierzchni nie lır 10-ciu 
w całej literaturze światowej. Swiadczy to pośrednio 
o dużej złożoności tego problemu i potrzebie kontynuowania 
badań w tym kierunku, 

Celem niniejszej pracy doktorskiej jest przeprowadzenie 
systematycznych badań nad rolą topografii powierzchni - 
w adsorpcji na powierzchniach EE. pe Ado 
Dokładniej, badania te dotyczyć beds “eve Gê 26) 
powierzchni na proces adsorpcji monowarstwowe ول‎ 
Zlokalizowanej, z bocznymi 619:007 w fazie 
adsorpcyjnej (rozdział 9 TM IV) oraz adsorpcji 
mobilnej (rozdział V). : 
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Il. MONOWARSTWOWA ADSORPCJA ZLOKALIZOWANA : 


BADANIA EFEKTOW  TOPOGRAFII POWIERZCHNI 
PRZY UŻYCIU PRZYBLIZENIA BRAGGA - WILLIAMSA 


hi.) Ogólne uwagi o modelach powierzchni heterogenicznych 
i zlokalizowanej adsorpcji monowarstwowe j czystych 


gazów na takich powierzchniach. 


Jako ilościową charakterystykę energetycznej 
heterogeniczności powierzchni przyjęto w teorii adsorpcji 
tzw. różniczkowy rozkład energii adsorpcji. Rozkład ten 
oznaczany przez X(E), definiowany jest jako OM/OE و‎ 
gdzie M jest ilością centrów adsorpcy jnych na 
powierzchni adsorbentu. Oczywiście funkcja X(E) 
charakteryzuje ilościowo jedynie globalną heterogeniczność 
powierzchni nie dając żadnych informacji, co do topo- 
graficznego rozkładu tej heterozeniczności, 

Rozważmy na początek model płatowy powierzchni. 
Powierzchnię adsorbentu rozpatruje się jako zbiór 
nieoddziaływujących, homogenicznych płatów w obrębie 
których istnieją identyczne warunki adsorpcji. 

Dla takiego modelu, fundamentalnym równaniem jest 
całkowa postać globalnej izotermy Gegner See 
V(p T) = m.” | (2.1) 
A 
zdzie V(p) jest izotermą adsorpcji globalnej, wyrażona; 
w jednostkach względnego pokrycia powierzchni, © Cp, Te) 
równaniem izotermy opisujycej adsorpcje lokalna na 


płacie powierzchni charakteryzuğacym się energia 205 jù 
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równą E, X(E) - znormalizowana do jedności funkcją 
rozkiadu energii adsorpcji, natomiast A dest 
DM DOR wartości energii adsorpcji istniejących 
realnie na powierzchni adsorbentu, W pracy będziemy 
przyjmowali, ze EE(E_,00) tje, że minimalna 
energia adsorpcji wynosi Bye 

Dla opisu układów adsorpcyjnych z powierzchniami 
heterogenicznymi, proponowano wiele analitycznych 
wyrażeń określających w przybliżeniu rozkład energii 
adsorpcji. Najstarszym z nich, jest chyba rozkład 
exponenc jalny/ Y , mający swe uzasadnienie w teorii 
tworzenia powierzchni رز‎ mee | 
Wojciechowski i współautorzy pokazali/ ^38 . Ze quasi- 
eksponencjalny rozkład Maxwella-Boltzmana opisuje 
heterogeniczność powierzchni bardzo dużej liczby 
układów adsorpcyjnych. Popularny rozkład gausowski, 
szczególnie przydatny w analizie danych adsorpcy jnych 
z obszaru niskich pokryć, wprowadzony został przez 
Oder Pröby vprovadzenia innych postaci funkc ji 
rozkzadw 37 nalezy uwazaé raczej za nieudane 
ze względu na brak ogölneğ n OE r 
Pomimo wielkiego sukcesu, obydwa: exponenc jalny 
i gausowski rozkład mają swoje słabe strony, na ktore 
, w wielu pracach zwracano niejednokrotnie uwage 13,147 
Należy podkreślić, że wszystkie te rozkłady znajduję 
swe racjonalne uzasadnienie w wyprowadzeniu z dańych 
eksperymentalnych przy zastosowaniu różnych technik 


matematycznych. Dla przykładu, jedną z metod bardzo 
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często używanych jest metoda transformacji Stieltjesa 
/15-17/. Jednakże brak zrozumiałego fizycznego 
pochodzenia wspomnianych funkcji rozkładu, mógł 

budzić wiele wątpliwości. Niedawna praca Cerofoliniego/ 19/ 
próbuje wyjaśnić ten problem. Jego hipoteza jest - . 
ulepszona wersja dobrze znanej w heterogenicznej 
katalizie, reguły "teta" 2f. Zgodnie z ty hipoteza, 
miejsca adsorpcyjne sa przypuszczalnie formowane jako 
rezultat cieplnego nieuporzadkowania. Stąd też, rozkład 
ich charakteryzowany jest temperaturą Te» w której 

ciało stałe jest preparowane. Centra adsorpcyjne 

uważano za defekty strukturalne wykazujące pewną "extra" 
energię U. Rozsądnym jest chyba twierdzenie, ze ich 
rozkład między "extra" energie U, charakteryzujący 
powierzchniową niejodnoroduość, jest w przybliżeńiu 
exponencjalny, jako prosta konsekwencja rozkładu 


Boltzmana 


M exp(-U/KT,) 3 
50 — — (2.2) 


W dawniejszych pracach/ 2, "extra" energia U była: 
identyfikowana z energiq adsorpcji E. Cerofolini, 
podchodząc bardziej realistycznie, rozważał U jako 
funkcję E. Zakładając, że funkcja U(E) i jej wyższe 


pochodne są ciągłe, można ją rozwinąć w szereg Taylora 


U(E) = 20,8 - ELE | (2.3) 


gdzie ER jest najmniejszą energia adsorpcji. 
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Pozostając tylko przy dwu pierwszych członach tego 
rozwinięcia i uwzględniając równanie (2.2) otrzymuje 


się następujące wyrażenie na funkcję rozkładu energii 


adsorpcji 
a,+28,(E-E_) a.(E-E_) a (E-E ye 
X(E) = (QM) Sh) 5 — ھ2‎ ET, o exp(— 2 一 和 o 
| (2,4) 
Jesli E-E 2? a,/a, wówczas 
2a (E-E) a.(E-B.)” 
X(E) = gy (a (2.5) 


Rozkład ten koresponduje z globalną izotermą adsorpcji 
. Dubinina-Raduszkiewicza (DR)/ 1% . opisującego wiele 


układów die کر شا‎ la al La 


przy niskich pokryciach‏ و 
powierzchni.‏ 


Jeśli natomiast E-E << SA مه‎ to otrzymuje się, 


CZ ) ظ‎ 


| tzn., rozkład korespondujucy z globalną izoterma 
| Freundlicha/ 19/, 

l W przypadku silnej heterogeniczności powierzchni, 

| szerokość funkcji rozkładu X(E) jest zwykle dużo 


| większa od kT, tak że z dobrym przybliżeniem możemy 


| całkowity proces adsorpcji rozważać jako proces 
skokowy, zachodzący w sekwencji male jących energii 
adsorpcji. Dlatego też, adsorpcja w obszarze niskich 
pokryć powinna być charakteryzowana wyższymi energiari 


adsorpcji, tj. przez warunek E-E 22 84/85, co prowadzi 
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Go opisu całkowitej adsorpcji równaniem DR . 


a 
Vpę(P) = exp(- gf, (eTinC/K*))*) (2.7) 
gdzie 
K’ = K exp(-E /KT) (2.8) 


zaś EK jest stałą z równania Langmuira O. (p,E), 
uzytego w charakterze izotermy lokalnej 


Q (aE) = (1 + Kexp(-B/kT)/p)"! SH (2,9) 


Jednakze, izoterma adsorpcji Dubinina-Raduszkiewicza 
Jest przy bardzo niskich pokryciach zbyt ogólna 

i nieczura na rzeczywisty charakter adsorpcji lokalnej 
z powodu chociażby tego, że w obszarze niskich pokryć 
powierzchni równanie Langmuira jest granicznym 
przypadkiem zarówno równania Fowlera~Guggenheima 

jak i Hilla-de Boera. 

Zauwazmy, ze topografia powierzchni nie vpiyva na proces 
adsorpcji, dopóki nie zaczna odgryva6 roli wzajemne 
oddziaływania pomiędzy zaadsorbovanymi molekułami. 
Oddziaływania te są efektywne w re jonie średnich. 
względnych pokryć powierzchni, kiedy to adsorpc ja 
scharakteryzowana jest przez warunek E-E << BA ane 
Odpowiednią funkcjy rozkładu energii adsorpcji jest 
wówczas z punktu widzenia teorii Cerofoliniego, | 
rozkład exponencjalny dany równaniem (2.6). -Rozkład 
ten, Łącznie z równaniem Langmuira jako izoterma 
lokalną, prowadzi do globalnej izotermy adsorpcji 

P i oba. 
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gdzie | 
r = SAIT, `` r€(0,1) (2.11) 


r - jest tzw. parametrem niejednorodności powierzchni., 


Bardziej ścisłe obliczenia wskazują, że równanie 


całkowe (2,12) 
8 r(E-E,) 4 
Ve(p) = | gy ود‎ pr) (U + Kexp(-B/x1)/p)” "az 
Ë, 
daje dla p<K’ 
Vp(p) = ا‎ + RTE ہے‎ (2.13) 
1-1 
natomiast dla p>) K?’ 
Vę(p) = 1 - 2 Ck 3 (2.14) 


i-1 
Równanie (2.10) jest dobrym przyblizeniem dla Vp(p), 
gdy szerokość funkcji rozkładu X(E) jest dużo większa 
od kT, który to warunek jeśt spełniony dla większości 
powierzchni heterogenicznych, 

Zanotujmy jeszcze odnośnie granic całkowania (E,,00), 
że taki wybór zakresu energii 050226 1 był dyskutowany 
przez Rudzińskiego i Nojeiechonskiege/ ER ; Wykazali oni, 
ze dolna granica całkowania Bo jest mocno zviazana 
ze stata Langmuira K, i niezalezne jej określenie 


z ciśnieniowych zależności danych adsorpcyjnych jest 
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niemożliwe. Hobson/ 27/ udowodnił natomiast, że ścisły 
wybór górnej granicy całkowania jest ważny jedynie 

przy bardzo małych pokryciach. Problem ten został także 
ilościowo przebadany przez Jarońca, Sokołowskiego 

i Rudzińskiego/ 2% , 

Provadzac zatem rozwazania w obszarze umiarkowanych 
pokryć powierzchni możemy bezpiecznie 6 (+00) 
Jako górną granicę energii adsorpcji, co jest także 
dogodne z punktu widzenia obliczeń matematycznych. 

W tym rozdziale zostaną przedstawione systematyczne 
studia nad wpływem topografii powierzchni na proces 
adsorpcji monowarstwowe j, zlokalizowanej z bocznymi 
oddziaływaniami między zaadsorbowanymi aR 
Scharakteryzujmy krétko model adsorpcji lokalnej przyjęty 
przez nas do badań wpływu topografii w tym rozdziale, ` 
Jest to prosty model Fowlera-Guggenheima, którego 
główne założenia możemy streścić następująco: 

a) na jednym miejscu adsorpcyjnym może być 

zlokalizowana co najwyżej jedna molekuła, 

b) każda ¿aa ds اد سا‎ molekuza oddz LAŁ je 

z pewną liczbą innych molekuł, które są 
zaadsorbowane na miejscach adsorpcy jnych 
najbliżej sąsiadujących, 

c) molekuły zaadsorbowane na miejscach najbliżej 

sąsiadujących z dana zaadsorbowang molekułą, 
są rozmieszczone na nich w sposób zupełnie 
| przypadkowy ° 
Przy tych założeniach, termodynamika statystyczna 


prowadzi do następującego równania izotermy adsorpcji 
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E+4kq O; 


© چم‎ (p E) = (1 + Kexp pp 2p) 


gdzie T. jest dvuvymiarova baie sd krytyczna 
adsorbatu. Jej wartość zwiyzaná jest z liczbq najblizej 
Sąsiadujycych miejsc adsorpcyjnych (c) i energią (w) 
wzajemnego oddziaływania dwóch molekuł zaadsorbowanych 


na. tych miejscach 


Szczegółowa postać stałej K dest znana ze statysty- 


cznego wyprowadzenia równania Langmuira 


3/2 
K . CERETI KT cu ) ری‎ 
h 0 


gdzie f jest funkcją rozdziału dla wewnętrznych 
stopni swobody molekuły adsorbatu nad 
powierzchniy 
f jest odpoviedniq funkc jy rozdziału dla 
zaadsorbowane j molekuły, która w przybliżeniu 
"oscylatora harmonicznego" wynosi qyqyq, 
gdzie d,s MT q, SA oscylacyjnymi funkcjami rozdziału 
dla wibracji w xierunku x, y; 2 odpowiednio. 

Pytanie, czy energia adsorpcji wpływa na wibracy jng 
tunke je rozdziału q,, ma wiele kontrowersyjnych 
odpowiedzi, Na przykład: teoretyczne rozważania Hilla 2% 
dają pierwiastkową zależność od D, dla częstości 
wibracji molekuły w warstwie zaadsorbowane 2. 

Drain i uğra ^. nie biorąc pod uwagę omawiane j 


zależności, otrzymali zgodność danych eksperymentälnych 


(2.15). 
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i teoretycznych. Appe2/ ^ V wyjaśnił temperaturowa 
zależność danych doświadczalnych dla pewnych układów, 
potwierdzoną przez obserwacje Halseya/ 3% , przyjmując 
exponencjalną zależność q; od E, Hoory i Prausnitz/ ^^/ 
wykazali, ze zwi&zek a E nie vpiyva praktycznie 

na proces adsorpc ji poğedyüczych gazów, Steele/ >*/ zaf, 
był odmiennego zdania, krytykując pracę Rossa i Oliviera 
N/ 

N prowadzonych przez nas rozvazaniach, bedziemy 


przyjmowali, ze zmiany E nie vpiyvağu na 6 f 


11.2 Adsorpcja na idealnie płatowych. powierzchniach 
heterogenicznych | 


Idealnie płatową powierzchnią heterogeniczna, 
będziemy nazywali powierzchnię adsorbentu, na której 
rzeczywiście można wyróżnić jednolite płaty. W zasięgu 
zaadsorbowane j molekuły znajdują się zatem miejsca 
adsorpcyjne scharakteryzowane ty samy energią adsorpcji, 

W całym rozdziale drugim przyjmujemy model adsorpcji 
lokalnej Fowlera-Guggenheima, oraz exponencjalny rozkład 
energii adsorpcji. 

Przy założeniu idealnie pzatowej topozrafii powierzchni, 
całkowe równanie dla globalnej izotermy. adsorpcji ma 


następujący postać : 


R AkT r 
E ۱ 
p FG - 1% 
V^ (p) x E exp(- zm )1 + K'exp( سس‎ )/ p) 4E 
O 


j aci a vielko6 si siq 
gdzie wskaźnik (P) oznacza, że dana wielkość odnosi się 


do modelu płatowego. Indeks (FG) w dalszej części 


RAGE" 


| rozdziału będzie opuszczony. 


Całkując równanie (2.19) przez części, otrzymuje się 


ES" 

| VG) = VCD) + vət) (2.20) 
gdzie 
" Bkt O” 4-1 
VA = [1 * K'exp( — ep Al, (2.21) 

00 
V5(p) = 2 OF (2.22) 
0 
Podstawiajyc 
E + net. ) E) ım. E 225) 
gdzie 
Ë 
t, - Ə (Bİ, (2.24) 
=O 
1 wprowadzajuc oznaczenie 
K'exp(-t /kT)/p = exp(s/kT) (2.25) 
š E QF ! | staci 
możemy zapisać w następującej pos 
۱ 105) = (14 exp SDE (2.26) 
Carka v5(s) po wprowadzeniu nowej zmiennej t 
Przyjmuje postać całki Fermiego-Diraca - 
ə S-t | 
| | exp( Sa? 
Vó(s) = yr | 505,9) —— (2.27) 
(1+exp(7z")) 
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a o1. 
gdzie 
H(t,s) = exp(- grtet -akr @*)) | (2.28) 


Rozwijajac funkcje H(t) w szereg Taylora w otoczeniu 
punktu tas 
oc 


Li 52 
H(t) = LES (SB) (2.29) 
ہے‎ il dt tas 


Wprowadzając novu zmienną x=(t-s)/kT otrzymujemy 
00 
Za A UE i 
Va(s) 5 ہس‎ Irae 
| i! dt tes )(1:e” )(14e ^) 
-8 
kT 


(2.30) 
ke 


Dla carek wystepujacych w równaniu (2.30) przy jmiemy 
pewne uproszczenie. Stosownie do równań (2.24) i (2.25) 


mamy, że dolna granica całkowania 


Es. O (2.31) 
KT e E=0 


Wpływ topografii powierzchni będzie niewątpliwie 
najbardziej interesujący w obszarze średnich i wyższych 
pokryć powierzchni, kiedy to oddziaływania między 
zaadsorbowanymi molekułami odgrywają największą role. 
Przyjmując zatem w tym rejonie całkowitych pokryć 
powierzchni, że miejsca o najniższej energii adsorpcji 
Sà zapełnione w połowie, otrzymujemy zgodnie z równaniem 


(2.31) olen EE 

p F 
ds AS (2.32) 
(1+e*)(1+e7 X) 1 
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gdzie 


| 00 
1 8 (189) 
j=0 
lub wyrażone przez funkeje zeta Riemanna $(1) 
A, = Cheti) bi) (2.54) 


Np., A,=1/2,1n2,9%°/12,... odpowiednio dla n-0,1,2,... 
Ograniczajac się tylko do pierwszych trzech członów 


rozwinięcia (2.29) otrzymujemy (2.35) 


25.2 ظ 

HD )‏ پر ارہ چک اس رش ) = Vals)‏ 
‚dt =8‏ 

W szczególności, 


H(s) - (p/K*)Fexp(2rT /T) 


H'(s) = (p/K')'exp(2rT /T)(r/kT) (T /T-1) (2.56) 
Bees) = (p/K* JFexp(2rt /2) (z/kT)* (T /T-4)° 


gdzie przez H? i Hit oznaczono pierwszą i drugą 

pochodne funkcji H(t). 

ięcząc równania (2.21) i (2.35) wg. (2.20) i uwzględnia- 
Jac (2.36) dostajemy jako kańcowy rezultat, następujące 
równanie dla izotermy adsorpcji globalnej, na płatowo 


heterogenicznych powierzchniach 
(p = @P(E_) + 0.5(p/K')Texp(2rT./T)(1 + 2r1n2 


2 ۳ 
“(1,/7-7) + dE XA. m3) (2,37) 


一‏ 22 ہے 


11.5 Adsorpcja na heterogenicznych powierzchniach 
wykazujących topografię typu "random" 

Ten model topografii powierzchni heterogenicznej, 
zakłada że nie ma żadnej przestrzennej korelacji 
między centrami adsorpcyjnymi o tej samej energii 
adsorpcji. Rozkład miejsc adsorpcyjnych na dowolnym 
Wycinku powierzchni a więc i w obszarze oddziaływań 
molekuły zaadsorbowanej jest dokładnie taki San’ 

1 określony przez globalną heterogenicznoéé powierzchni 
(funkcję rozkładu energii adsorpcji X(E) ). 

Przyjmując, że potencjał średniej siły 
działający na zaadsorbowany molekułę ma postud/ "%, 
ew (p) (2.38) ˆ 
Czyli nie zależy od tego na jakim centrum adsorpcy jnym 
została ona zlokalizowana, równanie izotermy lokalnej 
(FG) można zapisać w postaci | 


E+4kP_V*(p) 


©" (p,E) = (1 + K’exp ——crum ypr! (2.39) 


Indeks (R) oznacza tutaj topografię "random". 


Równanie catkowe izotermy globalnej ma wiec postaé 


oo 
"toi = J Em exp(-rE/kT) C (y ,E)aE ECKE 
8 


Catkujuc przez części otrzymuje Się 


(p) = V4(p) + y a (p) کیٹ‎ (2.41) 


"əə 


gdzie 
vite) = OPE) : (2.42) - 
= s-E 
V5(p) = E TEN AD a Er dE (2.43) 
(1+exp (sar) ( 
0 


powyzej exp(s/kT) jest obecnie równe K'exp(-A4T V /T)/ p. 
Całka (2.43) może być rozviazana analogicznie 


do równania (2.27). postępując tak, otrzymujemy 


H(s) = (p/K*)lexp(rT V/T) 


li 


H (s) (p/K')'exp(arp V"/T)(-2/kT) (2.44) 


H"(s) = (p/K')'exp(arp V /T)(r/ kT)" 


Podstawiajye (2.44) do równania (2.35) a następnie. 
łącząc z (2.42) według wzoru (2.41) otrzymujemy 
postać izotermy adsorpcji globalnej na powierzchniach 


typu "random" 


VR(p) = ۵۶ (p,E ) + 0.5(p/K*) exp سی‎ VI /T) 


Z : 
(1 - 2rln2 + İr”). (2.45) 


II.4 Adsorpcja. na hetercgenicanych powierzchniach 


0 pośrednim rozkładzie toposraficznym mie,sc 
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Jest mato pravdopodobne aby rzeczywiste 
powierzchnie heterogeniczne wykazywały zupełną 
korelację miejsc scharakteryzowanych tą samą 
energią adsorpcji lub też jej zupełny brak, 

Omówione do tej pory modele topografii powierzchni: 
model płatowy i model "random" należy uważać jako 
przypadki skrajne. Uzupełnimy je w tym rozdziale 

tzw. modelem pośrednim, 

W celu ilościowej oceny wpływu ograniczonej 
(pośredniej) przestrzennej korelacji miejsc 
adsorpcyjnych, wprowadzamy tzw. lokalny, różniczkowy 
rozkład energii adsorpcji X4 (E)E), który def inl kieny | 
nastepujaco: 

Rozwazmy miejsce adsorpcyjne scharakteryzowane 
energią adsorpcji E. X,(E,,E) definiujęmy jako 
różniczkowy rozkład energii adsorpcji Ey» Ç 
na heterogenicznym piacie-najblizszym sysiedztwie 
centrum adsorpcyjnego o energii E. 

W przypadku idealnie pzatowego, topograficznego 
rozkładu, X (E, ,E) jest funkcğa delta Diraca 


X (E,,E) = ô (E,-E) (2.46) 


W przypadku całkowicie przypadkowego rozkładu miejsc -, 


adsorpcyjnych, funkcja X4 CE, ,E) jest niezalezna 
od E , czyli | 


Xi(E,,E) = X(Ej) (2.47) 
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Dia interesujacego nas przypadku pośredniego, 
proponujemy dla X,(E,,E) następujące wyrażenie 


exp( (Ei -E)/D) 


X, (5, E) = l 一 一 一 —U[s,.- 
(1«exp((E,-E)/D)) 


5 (2.48) . 
Jest to funkcja typu "dzwonowego" z maximum w punkcie 
Ej*E i szerekoáci D . Jakkolwiek D jest zależne 
od E و‎ będziemy tu rozvazali prosty przypadek, gdy 
ta zależność jest pomijana. Możemy także mövid, 
ze operujemy pewny średnią wartością D. 

Potencjał średniej siły działającej na zaadsorbo- 
vana molekułę od strony najbliższych, sąsiadujących 
z nią, zaadsorbowanych molekuł, jest równy 4KT V(p,E) > 
gdzie V(p,E) oznacza względne pokrycie najbliżej 
sąsiadujących miejsc adsorpcyjnych rozważanej przez 
nas molekuły (tj. względne pokrycie obszaru oddziaływań 
molekuły zaadsorbovane na miejscu określonym energią E. 
W najprostszym przybliżeniu V(p,E) może być 


reprezentowane przez funkcję 
00 


-E)/D) s 
E EM E Ahi: asawa dle ee (2.49) 
V(p,E) (Xp, 1 )مده‎ (E, LEDS 1 
0 


gdzie (2.50) 


Op,sg,) = (1 + Kexp((E «kr, OCo;E 2/ C-k1)))7 


C(E) = 1+exp(-E/D) jest współczynnikiem normalizacji 
lokalnej funkcji rozkiadu (2.48). 


Uzywajuc tego samego rozwiniecia dla całki Fermiego- 
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Diraca , mozemy wyznaczyé V(p,E) z pozadana 
dokładnością, W przypadku pośredniego rozkładu 
topograficznego wystarczająco dokładnym przybliżeniem : 
powinno być przyjęcie pierwszych trzech członów tego 


rozwinięcia. Zatem, 


SICH 
(GE) = ( Oly) + 2DIn2 EE + 


d ° Q(y,E ) : 
+ ¿MD ر چ‎ (2.51) 


Powyższe równanie zostało otrzymane przy założeniu, 
że C(E)=2. To uproszczenie zostało podyktowane 


oczywistym warunkiem zbieżności szeregu (2.51) 


lim V(p,E) = ©)2,8( (2.52) 
DO 


Równanie całkowitej izotermy adsorpcji v“(p), 


przybiera obecnie następującą postać: 


(2.53) 


eo 
E44kT cV م)‎ E) 


vip) = | Erexp(-rE/kT)(1+K'exp(—pqp——)/P) dE 


Indeks (M) oznacza tu założenie o pośrednim modelu 
topograficznym (moderate topographical model). 
Całkując przez części i przeprowadzając zamianę 
zmiennych 

d 
E + 4k? V(p,E) = t + t, ۱ (2.54) ` 


gdzie 


t, = (60 (VCP LE) ( F (2.55) 
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otrzynujeny 
Vs) = vital + vis) my ۱ )2.56( 
Powyze j 
V (e) = 1/(d+exp(s/kT)) | (2.57) 
^ S-t 
M 1 exp) 
V_(8) = H(t,s) يقل‎ dt 2.58 
5 3 (1+exp( zk) (۴ Säit 
0 
H(t,s) = expl- gg(t«t,-4kT V(t,s))) (2.59) 


zag exp(s/kT)=K’exp(-t /KT)/p . Dalej, V jako 


funkcja zmiennej t ma postać 


V(t,s) < © + 2Dln2 ux A ee + 
(1-4 O(1- 9)) 


2 LE 
et” 2 Q@(1-Q)(1-29 ) 
6 (1-4(2 /T)0(1-O )) de 


gdzie O(t,s) jest zdefiniowane równaniem (2.26). 
Uzywajuc tej samej matematycznej metody do wyznaczenia 


və jak dla vo i və , otrzymuje sie 


H(s) = (P/K?) expl (2rT DEE 


2 2 
H'(s) = B(s) (z/KkT)(2,/7-1) (148 GL ET)" (T /m)) 


(4-9 2) t 12 
H"(s) = H(s)( /kT)*(T /T Pela OA E 
S) = s)(r T- + 一 一 -一 -一 一 
: 1261-7 /T) L 
ARIS 1. (2.61) 
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Jako ostateczny wynik otrzymujemy wyrazenie dla 
globalnej izotermy adsorpcji na powierzchniach 
heterogenicznych o pośrednim rozkładzie miejsc 


adsorpcy jnych 


vip) = Gg) + o. 5(p/k')"exp((1+ln2 ae ر‎ | 


“t/t 


at (D/KT)CT, /T) I 
کہ ہیں‎ +(D/KT)-(T —— 


ar” (D/ kT)? (T/T) 2 1 
^. 
1201-1 /T) | 


Ka,‏ پا 


+ ¿ETE لي‎ 1-1( 214 )2.62( 


W granicy D—0 , X, CE; ,E) staje sie Bake aa delta 
Diraca, 0(E,-E) » co oznacza و‎ ze nasza povierzchnia 
jest idealnie pzatowa. W tym przypadku równanie 

dla vip) poprawnie redukuje sie do równania na VE(p) 
opisującego. adsorpcję globalną na płatowo hetero- 
genicznych powierzchniach. 

Analizując równania otrzymane dla trzech różnych 
topograficznych modeli powierzchni, سا‎ NAL) że ۷ 


P 1 


granicy T ,/T-0 izotermy V و ۷ و‎ ۷ dane równaniami 


(2.57), (2.45), (2.62) yay re Gaja 6 do 


Vp) = (1 + و‎ + (y K?) sy Lei (2.63) 
> 1-0 


czyli do równania, które otrzymujemy przy założeniu. 
modelu Langmuira adsorpcji lokalnej. Możemy zatem | 
powiedzieć, ze gdy wzajemne oddziazywania między 
zaadsorbowanymi molekułami zanikają, topografia 


powierzchni nie wpływa juz na charakter adsorpcji. 
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11.5 Izosteryczne ciepło adsorpcji dla różnych 
modeli topografii powierzchni heterogenicznych 1 
Izosteryczne ciepło ad:orpcji jest definiowane 
równaniem 
Oe T(V_-V_) (C9 (2.64) 
st g a "aT V ۱ 


66 Qst jest izosterycznym ciepłem adsorpcji, 

ə i VA ۹ molovymi objetosciami odpoviednio fazy 
gazowej i warstwy zaadsorbowanej, p - cisnieniem | 
równowagowym gazowego adsorbatu, V - wzglednym 

pokryciem powierzchni adsorbentu. 

Zakładając idealno66 fazy gazowej i zaniedbując V, ə | 
można globalne ciepło izosteryczne dla pewnej wartości 


: 26/ 
pokrycia V zapisać równaniem Clausiusa-Clapeyrana >> 


Q.D = kr OD (2.65) 


Poniewaz całkowite pokrycie powierzchni V jest 


funkcją p i T و‎ dlatego przy stałym V 


(9 V(p,T)/ 8 (1/k1)), nos 7 


S) E ہج‎ (2.66) 
VU EIE (9 V(p,7)/ 8 (ln p))q 


W oparciu o powyzsze równanie mozna obliczyó 
izosteryczne ciepła adsorpcji używając do tego 

trzech równań izoterm adsorpcji vP, yə wr 
otrzymanych dla trzech róznych modeli topograficznych 
powierzchni. 


a) Izosteryczne ciepło adsorpcji na idealnie 


płatowej powierzchni heterogeniczne]j 
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Qa- OF 


1 Gë p, T) Cm AE Re ا‎ E 
| 1-4(T /T)@ )1- O) 


ررك من )هد + 


/ )))7/ دوک (21n2-1‏ رہہ سی engt‏ ’29/8 + )007-0 | 


| ^. 6©)1-©( 
1-401 /7) ©)1-©( 


+ r(vF-@)) i (2:67) 


gdzie 

O- C Co, E,) 

| Ë (2.68) 
| Q = 2.5kT + E, + 4kt © | 


b) Izosteryczne ciepło adsorpcji dla powierzchni 
| heterogeniczne j z topografia "random" 
| . . (2.69) 


; 9 (۷۲-۵ ( + 


af (p,n = [Gc © لور‎ + rkm((QR-ln( 


+ (P/K?) وہ‎ art, رد‎ (152-2 nêr) [9(1- 0) + (8-97 


gdzie 
O) = © (p, ) 

(2.70) 
Q = 2.5kT + E, + 4kT VF (p) 


c) Izosteryczne ciepło adsorpcji na powierzchni 
he terogeniczne j wykazujycej pośrednią 


topografię miejsc adsorpcyjnych 


Wprowadzając oznaczenia: 


L= 4(r/T7)9(1- Q) | | (2.71) 
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Tt (D/kT) (2, /7) 


ZĘ s 2.72 
12(1-7_/7) Bet 


równanie na ciepzo izosteryczne e możemy zapisać 


następująco (2.73) 


M eme M 
quM, EE 
se PT 1-L(0 V/90 ) ` 2 


© )+21n2* 
-© 


(D/ kT) (T A (D/KT)1n2 


— M_@) + ۰ ert AT EE 
ZEWIY )(V Séi )"exp((2r fT) 1-2 


+1))X21n2(1+S((T /T)“-(D/KT)(2+T / 2(-1( ۳ m )1-0 (۰ 


(A4S(D/kT)(A-SCD/kT) میم‎ ym) 1 EM ey" @)) 
ERC d 1-L(9 ۷/2 © ) 


gdzie 


O = Gl) ۱ (2.74) 


E 122022 ai ə 1 
Ee 1" AE ANES > + göt (D/ KT) T 


LI 
C28(1- © ( (4(07,/2)(9*-29 0-1) (1-29)*) (2.75) 


qu | ا‎ 7 E + تہ‎ (D nati 


“ny 
.. , ©U-6)(1-20) 
(1-2(1,/7)001-0)) + MD T Lem e ppt 
(2.76) 
W granicy T,/1”90 , QE. ` Zeit r zeA 
i jeżeli D30 otrzymujemy poprawnie nəə ۱ 
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II.6 Ilustracy jne obliczenia numeryczne do modelu 


Bragga-Williamsa 


Najpierw przedyskutujemy zachowanie sie 


P R 


globalnych izoterm adsorpcji: vt, m ates 


Wszystkie te izotermy mozemy przedstawió 


w następującej, ogólnej postaci 


/ 


)2.77( کر( + نهر( = vür,‏ 


gdzie i = (PRM) 6 gU jest ogölnie rzecz 
biorąc pewny funkcjy parametrów: r, T/T, ۰ 

Innymi słowami, każde z równać opisuğucych adsorpcję 
globalną jest liniowe kombinacją opowie عمس‎ izotermy 
lokalnej dla powierzchni homogenicznej z najbardziej 
prawdopodobny wartościy energii adsorpcji E, : 

i izotermy Freundlicha (p/K’)*. 

Drugi człon równania (2.77) nakiada pewne 
ograniczenia na równania izoterm vC) Ep), które 
wynikają z kilku przez nas przyjętych przyblizei. 

I tak, równania izoterm mogą być oczywiście 
interpretowane tylko w tym obszarze umiarkowanych 
pokryć powierzchni, w którym stosuje się izotermq 
Freundlicha. 

Chociaż wyniki obliczeń numerycznych będą prezentowane 
w całym zakresie globalnych pokryć powierzchni 
(przedział (0,1)), należy jednak pamiętać, że w rejonie 
bardzo małych i bardzo dużych pokryć dają przybliżone 


informacje, 


= 

7 3 éi 

FIR Tr 
= maż | 4 TR r 


22. 
" ə LE " 
E "yi پیجا‎ 
. Yo babe, 
Jc R 
JAN : 2 
z ya 
s x u (TR 
- D 
es 
Ye » 
E d 5 ^ 
r“ 
4 5 
e e uu. 
= w > " a 1 SÉ 
Har vd” 


T 


s -‏ ابص 
E - ` Ji. 3 A |‏ | 
X Tx .. Eno*a i Lay‏ 2 : 
m‏ ' 22 ويم d‏ 


T. 1 f 
“ur DE 


ç 
añ 
EI 


yu. و‎ ' 1 
Asi. om. Ke 
۱ ara 


š : 1 
qu — 0 1 5 PLI 
f s ^ => kla 


1 


- 35 ə 


Systematyczne badania nad wpływem topografii 
powierzchni na proces adsorpcji zaczniemy od porów- 
nania globalnych izoterm adsorpcji dla dwóch. ۱ 
skrajnych przypadków topografii powierzchni: piatowe j 
i w pełni przypadkowe j.-Rozpatrywaé będziemy tylko 
obszar temperatur superkrytycznych. 

Biorac 7 /#م‎ 5 WyZnaczono V”(p) i VF(p) dla 
trzech różnych parametrów niejednorodności: 

r= 0.01, 0.1, 0.3 . Wyniki przedstawiono na rysunku 4. 
Rysunek 2 umozliwia inne porównanie tych krzywych, 
Wyznaczonych dla stazej wartoóci r=0.1 i róznych 
wartości T/T (T /T- 0.1, 0.5, 0.75).. 

Ogólny wniosek jaki możemy wyciągnąć na podstawie 

tych dwóch rysunków można streścić następująco: 
powierzchnie 2 topozrafiy "random" miejsc adsorpcy jnycH 
zachowują większe podobieństwo do powierzchni homoge- 
nicznych w porównaniu z powierzchniami scharakteryzo- 
wanymi topografią płatową. Topografia płatowa uwypukla 
więc wpływ niejednorodności energetycznej powierzchni. 
Przy małych pokryciach adsorpcja na powierzchniach 

z topografia "random" jest mniejsza niż dla modelu 
płatowego, odwrotna relacja zachodzi dla wyższych 
globalnych pokryć powierzchni. Wyjaśnienie tego może 
być następujące: 

Dla wysokich wartości energii adsorpcji E , mamy 

w ogólności C (tp), Ep) . A viec, przy małych 
ciánieniach adsorbatu gdy proces adsorpcji zachodzi 


przede wszystkim na miejscach wysoce energetycznych 
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adsorpcja na powierzchniach z topografia "random" 
będzie mniejsza niz na powierzchniach z rozkładem 
topograficznym piatowym. Efekt przeciwny będzie 
zachodził przy pokryciach wysokich, kiedy to adsorpcja 
przebiega głównie na miejscach energetycznie as ron 
i wówczas Q(p)qv*(p) ۰ 

Przejdziemy teraz do pośredniej przestrzennej korelacji 
miejsc o tej samej energii adsorpcji. 

Rysunek 3 pokazuje jak zmienia się kształt globalnej 
izotermy adsorpcji V(p) kiedy przestrzenna korelacja 
miejsc adsorpcyjnych zmienia sie od płatowej (D=0) 

do "random" (oznaczonej na rysunku symbolem 00 ( ۰ 
Oczywiście w tym ostatnim przypadku YE p) obliczone 
zostato z równania (2.45) . Wyniki przedstawione na 
rysunku 3 sa otrzymane dla r=0.01 i T /T-0.5 : 
Wartość r-0.01 została tutaj przyjęta aby bardziej 
uwidocznić konsekwencje założenia różnych modeli 
topografii powierzchni, które sy uwypuklane dla 
powierzchni silnie heterogenicznych. Na rysunku tym 
mamy logiczną kolejność krzywych związanych z różnymi 
modelami powierzchni. Jest ona oczywiócie odbiciem 
zmian przebiegu izoterm globalnych. Kolejność ta 

może być też prostym sprawdzianem poprawności 
teoretycznych rozważań dotyczących topografii po- 
wierzchni. Rysunek 3 pokazuje również skutki 
wynikające z ucięcia szeregu (2.51) na trzecim 
członie - wartości D/kT wyższe niż O.^ nie 


mogą być brane do obliczeń ponieważ przy wyższych 
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ciśnieniach w tym przypadku mamy vp) VE (p) š 
Oczywiście, najwyższe możliwe do przyjęcia wartości 
D/kT zależą od parametrów r i T ui . To ogra- 
niczenie mozna zmniejszyé biorac pod uwage odpowiednio. 
większy ilość członów rozwinięcia (2.51) + 

Rysunek 4 pokazuje jak globalna niejednorodność 
powierzchni wpływa na przebieg zlobalnych izoterm 
adsorpcji w przypadku pośredniej przestrzenne j 
korelacji miejsc adsorpcyjnych. Następny rysunek 
dostarcza podobnych informacji dla zmieniających 

się wartości (me : 

Rysunki 6-10 obrazują wpływ topografii powierzchni 
na izosteryczne ciepła adsorpcji. Wszystkie krzywe 
obliczone zostały dla temperatury T=100 K . Przyczyną 
dla której zajmujemy się badaniem ciepeł adsorpcji 
jest niewątpliwie fakt, że dostarczają nam one 
bardziej obszernej charakterystyki adsorpcyjnej 

w porównaniu z izotermami adsorpcji. Tak jak | 

w większości publikacji و‎ prezentujemy ciepła adsorpcji 
jako funkcje względnego pokrycia powierzchni. 

Rysunek 6 daje porównanie izosterycznych ciepeł 

dla dwóch skrajnych modeli topografii powierzchni 
płatowego i "random" .Wykreślone tu krzywe są często 
otrzymywane eksperymentalnie i ich interpre tac ja 

może być następująca: 

Początkowe zmniejszanie się ciepeł izosterycznych 
odpowiada obszarowi pokryć w którym wpływ zmniejsza- 


jące j się wartości energii adsorpcji przeważa nad 
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wpływem wzrastającej energii oddziaływań między 


| 
| 
| zaadsorbowanymi molekułami. Z drugiej strony 


opadające gałęzie tych krzywych stają się szybko 
| liniowe zgodnie 2 równaniem Fowlera-Guggenheima 
| dla superkrytycznej adsorpcji na powierzchni homo- 
| genieznej. Innymi słowami, w tym obszarze pokryć 
| powierzchni wpływ niejednorodności powierzchni 
| szybko zanika. Latwo zauważyć, ze dla 0 
| modelu powierzchni zmniejszanie sie ciepeł adsorpcji 
|| zachodzi dużo szybciej. Pozostaje to w zgodności 
| z wcześniejszymi badaniami izoterm adsorpcji, 
ظ‎ gdyż gwartowniejszy spadek ciepeł adsorpcji sugeruje 
| większą niejednorodność powierzchni, | 
| Rysunek 7. pokazuje, że w pewnych przypadkach wykres 
izosterycznego ciepła adsorpcji może mieć tylko 
| opadajacy gałąź, co potwierdzają także dane 
13 eksperymentalne. 
| Rysunki 8-10 ilustruja zachovanie sie izosteryeznych 
ciepeł adsorpcjiw różnych sytuacjach fizycznych, 
w układach z pośrednią topografią miejsc adsorpcy jnych. 
Najbardziej interesujycym wśród nich jest rysunek 8 


| będący podstawową ilustracją wpływu topografii 


1 powierzchni na izosteryczne ciepła adsorpcji. 
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III. WYZNACZANIS ROZNICZKOWEGO ROZKLADU ENERGII 
ADSORPCJI PRZY ZAŁOŻENIU ROZNEJ TOPOGRAFII 
POWIERZCHNI  HETEROGENICZNEJ 

111.1 Ogólne rozwazania w oparciu o metode 


przybliżenia kondensacy ¡nego 


W niniejszym rozdziale pokażemy jak założenie 

o różnej topografii powierzchni heterogeniczne ل‎ 
prowadzi do różnych postaci wyznaczanych funkcji 
rozkładu energii adsorpcji. Użyjemy w tym celu metody 
przybliżenia kondensacyjnego, która dla przypadku 
topografii typu "random" została sformułowana przez 
autora tej pracy , w literaturze! V. | 

Istotą tej metody jest to, że jądro całkowego 
równania (2.1) - rzeczywista izoterma lokalna O (p,E) 
jest zastępowane przez skokową funkcję ©, (PE) ۱ 


(tzw. izoterme kondensacy jna) postaci 


o, p< p, CE,T) | 
‚E) = (3.1) 


którą można zapisać w równoważnej postaci 


x AN EC E (PT) 
E) = | (3.2) 
O. i * , E >B (pT) 


gdzie E (PT) jest funkcja odwrotną do p, CE ,T) 1 
Zastapienie (J(p,E) przez C) Cp, E) upraszcza 


problem rozwiazania catkowego röwnania Fredcholma I-go 


| 

| 
I 
| 
۱ 

| 

| 


MTI L et 


EN 


rodzaju (2.1), sprowadzajac go do równania 
00 
V(p) = | X(E) 00 (3.3) 
| E (p T) 


Rozwiązanie jest następujące Y 


X(E) = -|-> v? (5.4) 


p=p, (5) 


analityczna postać p, (E), zależy oczywiście od 


| założonego modelu adsorpcji lokalnej. wyznaczana jest 


Wr ware A 


© (p,E) = 0.5 (3.5) 


ktöry to warunek otrzymuje sie na bazie prostych 
rozważań wariacyjnych dla funkcji jednoznacznych 
x o wartoóciach z przedziazu (0,1). 
4 Przyjmujyc model Fowlera-Guggenheima adsorpcji lokalnej, 
| równania izoterm lokalnych OP), OX(p,E>) łącznie 
| z równaniem (5.5) dostarczajy nam następujących relacji 


między ciśnieniem kondensacji p, 4 energia adsorpcji E 


| EE) = exp((B-E0+2R0,)/-RT) (3.6) 


P(E) = exp((E-E0+4RT_V*(p))/-RT) EY? 
korespondujacych 2 piatowym modelem topograficznym 
oraz modelem typu "random" . Powyzej przez E? oznaczono 


wielkość RT mK. 


1 


ei! 
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Se 


= UB: 


Wyrażając eksperymentalną funkcję V(p) jako funkcję 
energii przy pomocy równań (5.6) i (3.7) a następnie 
różniczkując zgodnie z równaniem (3.4), dostaje się 
dwie funkcje rozkładu XE(E) i xn (E) korespondujace 
odpowiednio z pratowym i beziadnym topograficznym 
rozkładem. W przypadku rozpatrywanego tu przybliżenia 
Bragga-Williamsa otrzymujemy następujący ogólny 


związek między tymi funkcjami 
۱ 93 
X ظ)‎ (E رر‎ 


at PZA (2:87 


X^ (BR) 2 - 
1:4RT X (E (E )) 


Z powyższej zaleznoéci łatwo można wydedukować, ze 
funkc ja rozkładu Dec? otrzymana przy założeniu ۱ 
modelu "random" powierzchni, bedzie bardziej spiasz- 
czona w porównaniu do Xt (E) و‎ pokazujyc tym samym 
większy heterozeniczno$é powierzchni. Y granicy T=>0, 
przybliżone fubkode rozkładu energii adsorpcji 
wyliczone w oparciu o równanie (3.4) stają się 
identyczne z dokładnymi funkcjami rozkładu efor EP RTR 
zatem efektywność metody przybliżenia kondensacy jnego 
nie może poddawać w wytpliwość podstawowej konkluzji, 
że brak przestrzennej korelacji miejsc adsorpcyjnych 
o tej samej energii adsorpcji zmniejsza wpływ 
niejednorodności powierzchni. Ten wniosek można także 
wycięgnąć z badań modelowych sparmuay 7 


Rozważmy jeszcze relacje odwrotne 


EP = -RTInp + E? - RT, (5.9) 


E. -RTlnp + D e ۸۳۲ )( (3.10) ` 
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Z powyższych równań otrzymujemy 


E" < EP + ZRT_(1 - 2N(p)) (5.11) 


W granicy 0ه لا‎ , metoda CA (Condensation Approximation) 
oszacowuje najwyższą wartość energii adsorpcji 
istniejącą na badanej powierzchni. Równanie (3.11) daje 
nam następujący związek między maximum energii adsorpcji 


EE obliczonej przy założeniu pratowego topograficznego 


rozkładu i maximum energii adsorpcji E az obliczonej 
. 
przy założeniu topografii "random": 
R pi P ١ 
Bet = Bast ORT, ER. 


Badając granicę V—*1, dochodzimy w sposób analogiczny 


do następującej relacji między minimalną energią 


adsorpcji ER, a 0 a 
P 
Cas = Emin م6‎ CHE) 


Tak więc, w przypadku powierzchni heterogeniczne) 

o bezładnym rozkładzie miejsc adsorpcy jnych, otrzymujemy 
większy zakres energii adsorpcji istniejących 

na powierzchni. Innymi słowami, wyznaczona funkcja 
rozkładu energii x" (E) jest nie tylko spłaszczona 

w porównaniu z Get ale także rozszerzona o wartość 
ART... Powyższe konkluzje można podsumować nastgpująco: ! 
Założenie rozkładu topograficznego miejsc adsorpcyjnych 
typu "random" wskazuje na większą dyspersję energii 
adsorpcji na danej powierzchni heterogeniczne j 

v porövnaniu z topografia ۰ 

Rezultat ten nie powinien nas dziwić. Wyższe wartości E, : 


będą ogólnie dawać O(p,E)>V(p). W konsekwencji, 


可 


if 
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2:51 ۱ | IDA 
przy małych ciśnieniach kiedy adsorpcja przebiega 
głównie na miejscach wysokoenergetycznych, adsorpcja. 
na powierzchniach z topografig "random" jest One jez 
niż na powierzchni z topografią płatową (przy tej 
samej globalnej funkcji rozkładu X(E) ). Wyższe 
ciśnienia adsorbatu będą dawać efekt przeciwny. 

Jest to dobrze znany fakt, że globalna heterogeniczność 
powierzchni podwyższa adsorpcję w obszarze małych 
ciśnień i obniża ją w obszarze wysokich ciśnień 
adsorbatu. Możemy zatem powiedzieć, że przypadkowość 

w rozkładzie miejsc adsorpcyjnych na powierzchni 
niweluje w pewnym stopniu wpływ globalnej niejedno- 
rodnoéci. Do otrzymania tej samej krzywizny izotermy 
całkowitej adsorpcji V(p) musimy w przypadku 
powierzchni z topografia "random" przyjąć istnienie 
większej globalnej niejednorodności. Innymi słowami, 
przy tej samej krzywiźnie izotermy globalnej, każda 
metoda będzie dawać większą globalną heterogeniczność I 
powierzchni przy założeniu rozkładu "random" miejsc x 
adsorpcyjnych. 

Metoda CA wykazuje kilka zasadniczo róznych cech 
w zastosowaniu do przypadku płatowego i "random" 
topograficznego rozkładu. Zależność (3.9) można 


zapisać w następującej bardziej ogólnej postaci: 
x. Las RTlnp + const. , (34,14) 


gdzie stara zalezy od modelu adsorpcji lokalnej "Ara 


i wpływa jedynie na przesunięcie x (m) wzgledem 


Tw 


! 
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energetycznej skali gdy przechodzimy od jednej 
do drugiej izotermy lokalnej. Podobnie zależność (3.10) 
można zapisać w bardziej ogólnej postaci: 


RTlnp + £(VCp( ER )) ) (3.15)‏ يدم 


w ktorej funkcja f(V(p)) zalezy od izotermy lokalnej 
i wywołuje zmiany w kstałcie kondensacy jnego rozkładu, 
przy przejściu od jednej do innej izotermy lok-lne]j. 
W przypadku idealnej fazy zaadsorbowane ول‎ kiedy nie 
istnieją oddziaływania między zaadsorbovanymi .molekułami 
topografia powierzchni nie wpływa na dane adsorpcyjne. 
Tak jest na przykład, gdy użyjemy izotermy lokalne j 
Langmuira, do której redukuje się równanie Fowlera- 
Guggenheima w granicy T_ 0. Wówczas xe (E) i ELE) 
powinny redukować się do tej samej funkcji rozkładu 
energii adsorpcji. Równania (3.8) i (3.11) pokazują, 
ze tak w istocie jest. 

W celu ilustracji obliczyliśmy funkcje rozkładu 
energii CE) 1 XA (E) używając dwóch popularnych 
równań empirycznych opisujących adsorpcję globalną: 


równania Dubinina-Raduszkiewicza (DR) /porównaj rów.2.'7/ 
Vpg (P) = exp(-B(PT1n( p/p, JI ۱ ۱ (5.16) 
i izotermy Freundlicha (F) يسن دي‎ dA M rów.2.10/ 

V,.(p) = Ap” | SATA 


gdzie B, Pp» A i c są parametrami charakteryzującymi. 
dany układ adsorpcyjny ٠ 


Przyjmując płatowy cozidal tocozraficzny otrzymujemy 
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nastepuğace wyrazenia na funkcje rozkzadu energii 


d um E 
— ——— a 5 


| adsorpc ji: 
| 
P 70707 
Xop (E) = 2B(E-E5")exp( -B(E-EDR)2) (3.18) 
d 
| p = (c/RT)exp(-c/RT)exp(E-B ) (3.19) 
l gdzie 
¿DR xo ١ 
EG we BT. ORT, - RTln(p_) (2.20) 
EQ = E? - 2RT مضود‎ 47°) (5.21) 


Funkcje rozkładu Ziel?) i XB(E) mogą być łatwo 
wyznaczone z równania (5.8) . 

Jako typowy, rzeczywisty układ adsorpcyjny, do którego 
stosuje sig równanie DR., wybraliśmy układ  azot- 

örafit w T=78 K, zbadany przez Kiseleva i Isirikyana/ 9/ ° 
Parametry DR dia tego ukiadu zostaxy opublikovane” Y 


6 cal.”mol”, 9×ط‎ mops 


i wynoszą: B=6.93+107 
Wartość 21:050 100 Tor została obliczona w oparciu 

O WZÓT Raba الله‎ ə Dla cvuvymiaroveğQ temperatury 
Krytycznej uzyto wartosci 1,765 K و‎ otrzymaną przy 
założeniu, ze zaadsorbovana warstwa molekuł azotu 
zachowuje się jak dwuwymiarowy gaz Van der Vaalsa. 
Naturalnie, ostatnie założenie jest zasadniczo kontro- 
wersyjne w odniesieniu do modelu Fowlera-Guzzenheima 
1 dlatego akceptacja przez nas = in om ma tylko 
i wyłącznie charakter ilustracyjny. Chcemy jedynie 
posłużyć się bardziej realistyczną wartością Ioe 


Używając tych parametrów  wyznaczyliámy XE (E) 3 Xp CE) 
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korzystajac z rövnaü (5.8), (3.17). Porównania tych 
krzywych mozna dokonaó na rysunku 5.1 . Z porównania 
tego mozemy wyciagnaó wniosek, ze topografia bowe Aiki 
jest bardzo ważnym czynnikiem determinującym własności 
adsorpcy jne powierzchni heterogenicznej. 
| Rozważmy obecnie równanie Freundlicha. Jako typowy 
d układ adsorpcyjny opisywany ta izotermy, wybraliśmy 
układ argon-jodek cezu w T-78 g/ 10/ | Parametry 
dla tego układu zostały również opublikowane” V 
i wynoszą: A=0.152 , c=0.623 . Wartość K obliczona 
według formuły Hobsona wynosi dla argonu 2.6 408 Por; 
Dwuwymiarowa temperatura krytyczna dla argonu, została 
ustalona w ten sam sposób jak poprzednio Pe K = 
Uzyvağac powyzszych parametröv wyznaczylismy : 5و‎ 
1 Git . Wykresy tych funkcji rozkładu pokazane sa 
na rysunku 3.2 . Na podstawie tego rysunku, de e 
musimy dojść do wniosku o dużym znaczeniu topografii 
powierzchni. Zaskakujący więc jest fakt, że tak ważny 
z punktu widzenia fizycznego czynnik nie był uważnie 


studiowany w przeszłości, 


III.2 Wyznaczanie różniczkowego rozkładu energii 
adsorpcji przy zastosowaniu asymptotycznie 


słusznego przybliżenia kondensacy jnego 


Rozwiązanie równania całkowego (2.1) ze 
względu na funkcję rozkładu X(E) przy założonej 


izotermie adsorpcji globalnej V(p) jest także 


Pis X 


> 
LI 
- - 
WOTIS 
ST 
= Ww... 3 3 f 
qı — y — | 
5 
RE 
J m aş 
te 


DEERE 92008 


tsniedo ه”‎ Zog": SLL. 


unu 


proste gdy izotermq lokalna (9) p,E) zastąpimy 


tzw. izotermą "asymptotycznie słuszną" ©) €p,8), ۱ 


E:  (p/K)exp(£/RT) p > p, CE,T) 
و‎ ١و‎ “m, j 


(3.22) 
۹ b p > p. (E Mu 


a więc kombinację izotermy Henriego i funkcji skokowej. 
Parametr p, jest jak poprzednio ciśnieniem przy 
którym zachodzi objętościowa kondensacja, pokrywająca 
całkowicie monowarstwą rozpatrywane miejsca adsorpcy jne. 
Nazwa "asymptotycznie słuszna izoterma" wzięła sie 


stąd, że. spełnia ona asymptotyczne warunki: 


1i un BD) += d 
lin Oye | (3.23) 
sd C) (p,E) = (r/K)exp(E/RT). 


Zastepujac O(p,E) przez ONG , równanie całkowe 
dla globalnej izotermy adsorpcji przekształci się 


do postaci: 


E oo 
C ` š 
V(p) = | X(E)(p/K)exp(E/RT ل(‎ 2 X(E)dE (5.24) 


0 Bo 


M 


a avo 
Rozwiązaniem tego równania jest mok. ar / 


Ki n 
p APA R A, A 
av QE “p=p (E) “OBS p=p.(Ë) 


gdzie związek między ciśnieniem a energi4 dany jest 
jak poprzednio tj. równaniami (5.6) i (5.7) odpowiednio 
dla modelu powierzchni płatowej i "random". 


Równanie (3.25) można również zapisaś w bardziej 


Ze 
P 
سے‎ 


l 
- 


+ 
, 
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ogólnej postaci, korzystajac ze wzoru (5.4) na rozkrad 
energii adsorpcji w zwykłym przybliżeniu kondensacy jnym, 
mianowicie 

(i),p(i) 


x CO G2) = xg) * E. niai - (3.26) 


gdzie indeks i-PvR و‎ odnosi się odpowiednio do 
powierzchni płatowej i "random". 
Dla topografii typu "random" otrzymuje sie, 


QXR (DxD. (AEF, og?) 


— W umur": d ا‎ 
"۵ g^ (1+4RT_XĘ(E" ))* 


Z równania (5.6) 1 (5.7) otrzymujemy dalej, 


> zl 1 


= U : (3 28) 
Oh) aT XE (E^) 


Ostatecznie wiec, 


(329) 
xPP) ی80‎ ۳/۵ EP) 
O AA حي ل‎ 
s VART XECE) ¿ (1+4RT_XQ(B )) 
RAR 
e Rat) 
- ar) + Ga EISE DT te y 


xc ay) 


Bazując na wynikach poprzedniego podrozdziału (równania 
(3.18), (3.19)) tj. na wyliczonych funkcjach rozkładu 
energii adsorpcji w przybliżeniu kondensacy jnym dla 
układów opisywanych izotermą globalny (DR) i (P) , 
ec sə I różnych modelach topograficznych 
powierzchni, wyznaczymy funkc je rozkładu energii 


adsorpcji w przybliżeniu kondensacyjnym asymptotycznie 


اب 


| j P - 57 m 
słusznym. ak, dla modelu płatowezo powierzchni 


otrzymujemy 


UJ 
— 
1 


2BRT(1-2B(E- BA?) *)e::p(-B(F 8 be + 2B(E E 5. 
-exp( -B(E-ED3)2) (DR) © (3.30) 


P 
EL 


(1-0)(c/RT)exp((-c/RT)(E-E)) . (€) (3.31) 


Funkcje rozkładu TAE) dla modelu powierzchni typu 
random" mogę być wyznaczone przy użyciu równania (3.29). 

E | Analizując równanie (5.26) można zaobserwować, 

x ze dla X (E)>O funkcja X4(E)» X. (E) , natomiast 

| dla X ,(8)4 O zachodzi relacja odwrotna tj. PROS < X (E). 

Należy również zauważyć, ze 'X (Œ) jest funkcją 
znormalizowans do jedności  jeóli korzysta się 


| ze znormalizowane ij postaci funkcji Xx (49 i spetniony 


jc neste pugacy "aura 

09 

| ox AXE. | | 

i = 一 Xo (Eo) zə) ER (3.32) 
pou | | 

l E 


Wyniki obliczeń numerycznych ilustrujęce wpływ. 
1 założonego modelu powierzchni heterogenicznej na rozkład 
energii adsorpcji w przybliżeniu asymptotycznie słusznym, 
są przedstawione na rysunku 5.5 . 
Funke je rozkładu xt (E) y XE) obliczono na podstawie 
| równań (3.30) i (5.29) , przyjmujęc parametry układu (DR) 
|  azot-grafit w Te78 K. Zzodnie z przewidywaniami otrzymuje 


` gig większą dyspersją wartości energii adsorpcji dla 


Kä? 


| 
} 
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powierzchni typu "random" w porównaniu z powierzchnis 
płatowy, przy załozeniu w obu przypadkach tej samej 
izotermy globalnej (DR) .W ten sposób uzyskano 
potwierdzenie wyników otrzymanych w poprzednim rozdziale. 
W przypadku funkcji rozkładu x (E) ,و‎ w Obszarze 
najmniejszych wartości energii adsorpcji pojawia się 
niesymetryczny pik czego nie obserwuje Się w zwykłym 
przybliżeniu kondensacyjnym. Podobny kształt funkcji 
rozkładu otrzymali Drain i Morrison/ 15/ z kaloryme- 
trycznych danych eksperymentulnych dlu ukłudów urgon- 
rutyl w T485 K i tlen-rutyl w "100 K. 

Na rysunku 5.5 wykreślono dla porównania z krzywymi 

X (E) , funkcje rozkładu energii adsorpcji X;(E) 

i x (5) uzyskane nu buzie zwykłego przyblizenia 


kondensucy nego. 


TOPOGRAFIA PLATOWA 
一 一 一 TOPOGRAFIA RANDOM 
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 Rys. 3.1 Funkcje rozkładu energii adsorpcji dla układu — 


azot - grafit w 78 K opisanego równaniem DR 
| 948 I x A 


przy założeniu topografii płutowej 1 "random" 


powierzchni heterozenicznej. "S$: ul, 
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TOPOGRAFIA PLATOWA 
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FUNKCJA ROZKŁADU ENERGII ADSORPCJI 


.0 2.0 
ENERGIA ADSORPCJI W KCAL/MOL 


Rys. 5.2 Funkcje rozkładu energii adsorpcji dla układu 
argon - jodek cezu w 78 K opisanego równaniem 
Freundlicha przy założeniu topografii pzatowej 


i "random" powierzchni heterozeniczne j. 
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IV. ZASTOSOWANIE PRZYBLIZENIA 4UAST- CHEMICZNEGO 
FOWLERA-GUGGENHEIMA DO BADANIA ADSORPCJI. 
NA POWIERZCHNI HETEROGENICZNEJ O 20230 
ROZKŁADZIE MIEJSC ADSORPCYJNYCH , 


IV.1 Ozólne wyrażenia dla funkcji termodynamicznych 
w superkrytycznym obszarze temperatur. 


Kanoniczna suma stanów dla układu N cząstek 


zaadsorbowanych na M miejscach adsorpcyjnych 


powierzchni heterogenicznej ma 8 


i 


2(N,M,T) = 2.2 (4.1) 
(Nj 


NT. NM ۱ 
/ | 11 | | N 
N, E 2. g(N44)exp( — ( [ [2 exp(E,/11)) 1 
N44 
gdzie: N; — liczba czystek zaadsorbovanych na M. 


Miejscach adsorpey jnych, scharakteryzowanych energią 


i. UE | Bà‏ ہد 


AQ _ liczbs konfiguracji z N,, parami 
żaadsorbowanych cząstek na miejscach adsorpcy jnych 
Najbliżej sasiadujycych, przy określonym rozkładzie 5 


IN} - In, NO AO 2 E (N44) TT) s 


i | 1 


ty molekularna funkcja rozdziału cząstki ` 


żaadsorbowanej na miejscu określonym energia adsorpcji Ej; - 
Do Obliczenia AEPD notrzebna jest dokładna 


£ 


informacia o sposobie rozmieszczenia centrów 


"A-a 


| 
1 


— uğ: 


adsorpcyjnych na powierzchni adsorbentu. 
Bi))/ V , rozwazajec powierzchnie heterogeniczna 
o bezładnym rozkładzie miejsc adsorpcyjnych, otrzymał 


ogólne wyrażenia na globalną izotermę adsorpcji V(p), 


energie wewnętrzny Y i energię swobodną Helmholtza A. 


Rozważania Hilla można streścić następująco: 

Liczba konfiguracji By Z DS. ال‎ 10828, 1) stanovi 
pewny, część liczby sposobów Bu = MI/C(NICM-N)!) و‎ 
rozkładu N czystek na M miejscach. Ponieważ nożkzAdł 
miejsc o różnej energii adsorpcji jest zupełnie 
przypadkowy, to liczbę By, konfiguracji można uważać 


za próbę losowy Z Boe WÓWCZAS 


&(344)/GCR44) F Pei Da, (4.2)) 


sdzie G(N44) jest liczba konfiguracji dla dane j 
wartości parametru N44: 
horzystając z powyższej zuleżaości 5(N44) od GCN44) 


można zupisać, ze 


(4.5) 
by 
N. 
[ =: > | -= 
SE N44 


Suma po prawej stronie jest take Sama jak dia 
powierzchni homəgeniczneğ ۳ ka herr SAR parametrów 

N, M, T, a nie zalezy od sposobu podziału N na N. e 
Zastyypienie tej sumy jej maxsymalnym CZŁONEM 

w punkcie NA A w przybiiżentu Guasi- chemicznym, 


4 D = > P 
b ss sna ödə” Yen OQ Let 
prowadzi. dg nas Span "یی‎ 


—a وس‎ — r 


Lé o 


B, 


N. N. 
zn, = G (NĄ, )exp(-wN2 4/ KT) SCH VG ار‎ km) tiy 
1 
("i | N, d 
۔‎ expo | (İLİ) E expen) (4.4) 
i ı 
gdzie 


(F-1+2V)(1-V 
1 ) in( 1" + 1-21 nl 


1 
YC = el «(kT/w)(V1n( , 
2 ۱ (F+1+2V ۷ (F+1)(1-V) J 


1/2 
F = (1-4(1-exp(-w/kT))(1-V)V) . (4.5) 


Dla znalezienia rozktadu Ny czustek pomiedzy M. 


miejsc adsorpcyjnych, dającego maksymalny Człon Ly 


“r 
którym następnie zastępuje się kanoniczną sume erte 
Z(N,M,T) , Hill zastosował metodę nieoznaczonych 
czynników Lagrange'a, "maksymalizujac" ln E wzgledem 
Ni przy ustalonym N . Otrzymał w ten sposób wyrażenie 
określające liczbę molekuł N, zaadsorbowanych ' 


na M. miejscach adsorpcy jnych 


M. 
1 


"dcum. تھا‎ c (4.6) 


i 1+exp(a-E,/kT)/f; 


gdzie (a) jest czynnikiem nieoznaczonym. 
Podstaviağuc równanie (4.6) do równania (4.4) 


otrzymuje się 


| (4.7) 
InZ(N,M,T) = aN-(wJt/kT)s رج‎ M, In( tat exp(S,/kP-a)) 


1 


m 
) + 


Różniczkowanie (3 1n2/ ON) daje potencjał chemiczny | 


عه سض ج و A‏ 


af rg 
——” 
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| Skradnika w fazie adsorpcyjnej Pa 


=- M /kT = a - cw/2kT - Mp iniu. Ads: Ur 4.8 
| Za 2 (F+1-2V)V M A 
| W stanie równowagi a7 Be (Hz - potencjał chemiczny 
| gazowego adsorbatu) . Zakiadajgc idealność fazy 

| objetoóciowej, mamy 


3⁄2 
- | KT = o... 3 70080 (4.9) 


gdzie p oznacza cignienie równowagowe adsorbatu 
| zaś K jest stałą Langmuira. 
Z równań (4.6), (4.7) i (4.9) otrzymuje sie 
Wyrazenie na izotermę „lokalną ©, = NM. , dającą sie 


zapisać w postaci równania typu Langmuira . 
3 i 
Q,;(p,T,E,) = (1 + K'exp( -E,/kT)/p) (4.40). 


gdzie 


exp(w/kT)(F-1+2V)(1-V) c/ 2 
K“ = K: (—— — TVW 1 (4.11) 
(F+1-2V ۷ 
Przy założeniu ciygłej zmiany wartości energii 
adsorpcji na powierzchni ciała stałego, wielkości 


termodynamiczne wyrażą się w formie całkowej przez 


następujące równania: 


۷ = | X(E)QO(E)de (4.12) 


A 


dla izotermy giobalnej Y, 


m 


= E ہے‎ EE _ 


~ 66 = 


1 A = Mk? | Juste Bi - av) + wj (4.13) 


dla energii swobodnej Helmhcltza A و‎ 


2 (v3 dlnf E ` 
E - MkT ¡ESTE EE OF + wot (4.14) 


aT 


e 


dla energii wewnętrznej E. 
Dla tak zdefiniowanych tu funkcji termodynamicznych, 
pojemność cieplna warstwy zaadsorbowane j C. 


, 


będzie mieć postać 


dlnf ， ہے‎ d لضا‎ ۱ 
Cy = ie | pca cama m kT. ep), + 
۵ 
dlnf, ام‎ 9-7 aX | 

- ACEI ET ع‎ 0... "Cam, (4.15) 

A 
gdzie 
— 2۷ (1-V) zą Se 
JC - O,5cN(1 = “ün | (4.16) | 
2ص‎ 
( . 2New 1-V) exp(-w/ Kü Déi 
“an”? kT F(F41)^ (4.17) 
Z równania (4.6) E 

d Ç E Lt dinf < 
) A - - © )1-9( (gg + ar, - (m) (4.18) 


Zależność (a) od V i T otrzymuje się z rownan 
(4.6) i (4.12) . Rözniczkowanie równania (4.12) 


względem T , przy stałym globalnym pokryciu ۷ و‎ 


= a " zi » 
= d | . 5 
تو‎ " zB سے‎ - m = ال سح‎ -— "LL. v رھ‎ >. ° 
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provadzi do: : . š 

| 7007 

ha ROGER, Zu CES 
A OT 0 "e d 


Stad, uwzgledniajac równanie (4.18) otrzymuje sie 


“Oa NOS z 9 (1- © )dE 
FR. ALLEE (4.20) 
=P któe 

A^ 


Aby wyprowadziá ogólne wyrazenie na bilne 
lzosteryczne Hat nalezy przede wszystkim uwzglednié 
fakt, ze izoterma lokalna @(E) jest funkcją para- 
metrów p i T także przez wielkość globalnego 
pokrycia V(p,T) , od której zalezy ezynnik nieoznaczony 


(a) . Sruszne sa zatem röwnosci: 


| 0 
SAD AOS (4.21) 
d T p Qv p,T oT p OT p,V 
| 0 © © 2 IV ٠ e O. - 2 
ہت‎ + a ) T.p p ölmp”n f “SIND my 0 pa 
oraz 
OV 
Os se np) ain + ont (4.23) 


Na podstawie ostatniego röwnania mamy 


(9 v/ dT). 
st Wye IN 
"T 


q 


d 


D 
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Korzystając teraz z równań: (4.12), (4.21) i (4,22)  - 


otrzymuje sie 


(-9.V 53 4۵ 
OT = s (4.25) 
AER, dz 
d Oy dE 
Cn = E ə (4.26) 


TO a سج‎ S), سے‎ 
— 


i ostatecznie | - 


a ES Ink | E cw(F-1*2v)] REES"‏ © م 
op Es E‏ ہر سم + | ST, = -O( O)‏ 


or Ja‏ س 


D | | ‘ b 4 = “ 


Ostatnie dwa wyrażenia zostały znalezione przy użyciu ` 


nun 
równań (4.10) i (4.11) . -. qi o qeu E Lir. həb 
Ogólne wyrazenie na globalne, izosteryczne ad 4 "ne 


adsorpcji dla powierzchni heterogenicznej typu "random 


jest analogiczne do wyrażenia otrzymanego przez ہے‎ 


A 


- mmm 
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| R.H.van Dongena i J.C.P. Broekhoffa/ Y przy założeniu 


płatowego modelu powierzchni 


ROCZE 
22 | 
dat ED ———.—. —r.h , (4.50) 


po $$, 


& 


gdzie lokalna izoterma adsorpcji © (p,T,E) 


nie zależy od globalnego pokrycia V , 


IV.2 Ogólne wyrazenia dla funkcji termodynamicznych 
przy zastosowaniu exponencjalnego rozkładu 
energii adsorpcji 


W dalszych rozważaniach zaniedbywać będziemy 
wpływ E i T na molekularną funkcję rozdziału ۰ 
zaadsorbowane j cząstki. Rozkład energii adsorpcji  X(E) 


aproksymować będziemy funkcją postaci: EN. 


X CE گور‎ sin(4tr) — A ود‎ (4.31) | 
TAT (exp((B-E,)/kT)-1) ə. 

Rozkrad taki wynika z ogólnych rozważań Cerofoliniego ۷ 

Jako graniczna postać ogólniejszego rozkładu, słuszna y 


W obszarze árednich i wyzszych pokry¢ powierzchni. 
Koresponduje on z uogólnioną izoterma Freundlicha, 
Przy założeniu równania Langmuira na lokalny izotermg 
adsérpe ji/ ^, 
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oraz funkcję rozkładu (4.51), otrzymuje sie przez 
Wycatkowanie równania (4.12) następującą izoterme 


globalną : 


NK MERE exp(a-E/kD)/$,) (4.32) a 
| D 

Stud czynnik nieoznaczony (a) jako funkcja V i T Em 
E y1/T Tak, 1 

AVE EL/ kö, )ما و كيدا رك‎ ۱ (4.33) RP 
se 


Uwzzlędnienie wyrażeń na potencjał chemiczny składnika du 
| Ma 
° M ۳ ` ۰ ۰ و‎ r = 4 2 > J 
w fazie adsorpcyjnej Ma w równaniu (4.32) p rn 


do izotermy globalnej postaci . TRE 


1/r 


| " 
(Cu MUA V 
P = Kennt äh TE 


PIDE. ao e 
Z8piszmy jeszcze raz równania na wielkosci termo- 


dynamiezrne w nieco innej postaci, ktora Jesti 


dogodnie jsza do dalszych przekształceń | a 
| é رکم جس‎ 


woe i. 
MKTC-Vis ى‎ KT d MET ^ Er. 


m 
On 


am. — 
—k 


wa xakə i 


ہے 


۔718, _ 


edlnK cw(F-1+2V 
GN? ALA — ' (3.38) 


00 
o, = uk] | مرو 26 فجي دق‎ + QM , 
E, 


o0 
( B-E 
(J) قمرقےت ری×‎ )1-0(۵8 )* 
E 


O 


Ce (4.39) 


co 
۱ x(5)Ö(1- Ö)dE 
E 


O 


Wykonując całkowanie przy zastosowaniu exponencjalnej 
funkcji rozkładu energii adsorpcji (4.32) otrzymuje się 
następujące analityczne wyrażenia dla funkcji termo- 
dynamicznych jako funkcje względnego pokrycia 


powierzchni V 


(4.40) 
CWV V 
A = MkT(-Vlnf_ - VE /KT + gp + Yin(——n + 
y Vr 
L 81 AE dae CY1n( F-1+2V)(1-V ) it ) 
سے 06ج‎ * 923v ao 9 * | xps 
0 
Bass ma kr A E dl ) - w(i + 
T deu | 7 T + spi سی‎ +F 
1 
x yVr | x lax ) (4.41) 
| x-v VT s > 
0 ۸ | 
ei 1 6 ) h y V r7 1 qlz im “lax 1 
Rat = 1 + kT -kT n a "ZZ: 


o 


AV "də -VKT _ cw F-1:2V) (4.42) 
"way VT) Mr ۱ 


7. 


Z SZ n i ١ '‏ موہ کا 
RK il y_y lr R dë‏ مہ m vi c 268”V (1-V)‏ 
PERO con‏ "(۶+1(۴)۱۲) ۶ ` 
1 
Int y Vr ) R Vel xP Tax ET sin(q r) i‏ ۷) ۰ 
T-V E ) x-V VE 1 4 TT T‏ 


0 


00 


ev 1/r 1/ 1n”(1:x dx | 
ray /T)“ | 20 
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Całki występujące we wzorach (4.40-4.43) wynoszy 


Odpowiednio 
ut/r 
ae لقره دہ کن مس رہ‎ j 
| ^u dY De eor M o (4,44) ہے‎ 
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IV 25 Adsorpcja w obszarze kondensacy jnym s 
uwagi ogólne 


ei Vöröd wielu zjawisk krytycznych jakie mogą mieó. 

۱ miejsce w układach "zależnych" gdzie istnieją 

oddziaływania między molekułami zaadsorbowanymi, 
najczęściej badana jest kondensacja dwuwymiarowa. d 
Kilka teoretycznych modeli faz ma nı EH 


przepowiada istnienie przejść fazowych pierwszego. 


rodzaju. Wśród nich najbardziej popularnymi są: 

model adsorpcji mobilnej Hilla-de Boera/ + j model 
adsorpcji zlbkalizowanej Fowlera-Guggenheim / 9f. 
Teorie te prowadzą do zadawalajacego jakościowo opisu 
wielu rzeczywistych układów adsorpcyjnych w pobliżu 

i poniżej temperatury krytycznej. Jednakże od bardzo 
dawna obserwowano również pewne dramatyczne rozbieżno- | 
Sci między wartościami doświadczalnymi a przewidywanymi | 
teoretycznie. Tak jest na przykiad w przypadku dwu- 
۔‎ Wymiarowej temperatury krytycznej T... 

0 Dla dwuwymiarowej , mobilnej warstwy adsorpcyjnej, 

mn 11/9 i de د‎ M przewidywali, ze dla kulistych, 
izotropowych molekuz, stosunek dwuwymiarowej tempera- 
tury krytycznej Toe do trójwymiarowe J pria Fl aa 
krytycznej T c gazu objętościowego powinien być A) | 
równy 0.5 . یمان کمن ےج‎ adoptujac teorię celową e 
Lennard-Jonesa-Devonshire dla przypadku dwuwymiarowego, ion" 
doszedł do wniosku, ze powyższy stosunek powinien być "ARE 


nieco wyższy, tj. rowny 0.53 


j mei LL" E بد‎ 55 R 
"E om "p "a “a 3 i ig 


MI. 
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m 


ہو وراد 


Tymaczasem, obserwowane dwuwymiarowe temperatury 


krytyczne sa znacznie niższe” *+8/ . 205153137 one 


wyznaczone przy najlepszym dopasowaniu teoretycznych: = 


wyrazen do danych eksperymentalnych/ 27 1¢/ lub bezpo- 
średnią obs PSS AA (WR W odniesieniu do pomiarów 
izoterm adsorpcji, istnieją قاذ‎ meni mine fakty 
wskazujące, że stosunek ten jest z reguły mniejszy 

od 0.4 , Dowodzą tego np. pomiary pojemności cieplnej 
warstwy zaadsorbowanej Ne, Ar, Kr na Graphonie 
wykonane przez Steele i Karla? 134 

Pierwsza próba wyjaśnienia tego faktu pochodzi od 

de Ban , Rozważał on możliwość odpychania dipoli 
zaadsorbowanych molekuł przez jony sieci nach Ronse 
Jednakże wielkość indukowanych dipoli, wyznaczona 
MezaTeżnie/ TIAM nie popiera tego tłumaczenia. 


/16,17/ 


Również trzeciorzędowe oddziaływania ‘nie mogą 


służyć jako zasadnicze wyjuśnienie, jak to zostało. 


2 kolei przyjęte przez de Boera i Broekho£t a^ 18/ 


Prócz tego, zgodnie 2 oszacovaniem Everetta 1% 
Wpływ trzeciorzędowych oddziaływań był znacznie 
przeceniany w dotychczasowych badaniach, 
Ross i Olivier Y byli piervszymi, ktörzy postulovali, 
Ze dwuwymiarowe temperatury krytyczne, mniejsze od 
wartości teoretycznych mogą być w pełni wy jaónione i 
przez niejednorodność powierzchni. 

Heterogeniczno66 powierzchni wpływa w Sposób 
bardzo skomplikowany na krytyczne zjawiska powierzch: 


niowe. I tak, w pewnych przypadkach adsorpcji 
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¿reficie o dobrze rozwiniętej Stran turze warstwowej 
(extrafoliated graphi te )/ 29-07 obserwuje się wiele 
pionowych skoków na izotermie acsorpe ji. Sa one 
przypisywane dwuwymiarowej kondensacji w pierwszej | 
i wyzszych warstwach zaadsorbowanych. Tomy i Duva1/ 25/ 
wykazali eksperymentalnie, że skoki te stają się 
bardziej wyraźne przy przejściu od naturalnego (nie- 
jednorodnego) grafitu przez Sterling MT (3000 G) 
do exfoliowanego grafitu. 

W rozdziale tym będziemy interesować się 
zlokalizowaną monovarstva uwzyledniajuca oddzi My “asır 
boczne, w temperaturach subkrytycznych. 

4 tego punktu widzenia najbardziej interesujyce, 
teoretyczne rezultaty były opublikowane przez Willa’ Y 


[22-29 


Gordona/ ^9 . Harris: van Donzena i Broekhoffa 


/30,31/ oraz przez Gerofoliniego/?€-39/ , 

Hill, zakłada jąc topografię "random" centrów adsorpcy j- 
nych i przyjmujyc w obliczeniach model quasi-chemiczny 
Fovlera-Gussenheima, doszed do wniasku, Ze w pewnych 
arunkac! ndensgació تا اه‎ DALE | JWUZUGOMLOWO 
niz jak zwykle jednostopniowo. nicsex ten, zostai 
także potwierdzony przez symulacje komputerowe Gordona 
/26/. który badał adsorpcję gazów na keterogenicznej, 
heksagonainej sieci z uwzględnieniem oddzieływań 
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Obecnie, przyjmując te same założenia co na początku 
rozdziału, spróbujemy dokonać opisu adsorpcji 
na heterogenicznych powierzchniach, w temperaturach 


Subkrytycznych. 


IV.4 Punkt krytyczny - przypadek powierzchni "random" 


Otrzymane vyrazenie na globalna izoterme adsorpcji 
VM (równanie 4.34) prowadzi do przegieé krzywe j lnp-V 
| FE dla odpowiednio dużych wartości -cv/kT , przepowiada 
a Wiec pierwszorzedowe przej$cie fazowe. 
Wartości parametrów krytycznych: Va (tj. wzgledne 
| pokrycie powierzchni w punkcie przegiecia izotermy 
| krytycznej) oraz B. (dwuwymiarowa temperatura krytyozna) 


2 Można wyliczyć 2 warunków: 


E A teg, د‎ = "inp | du 
2 CZY A 0 Czyż h | (4.47) 


Uzyvağac do tego celu równania izotermy (4,34) 


dochodzi sie po zmudnych rachunkach, do następującego 


układu równań 


ra 
a 
3.4 


öl c(4-F) 1 M | E 
(—) = Ce 5 = O (4.48) “x 
UV in TT کید‎ x ? 


۹/2 (4, 49) ۰ 
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do równania drugiego, otrzymuje sie równanie dogocne 


dla znajdownia wartości Ve و‎ Które są jego pierwia- 


stkami 
f, +f, =0 | ` (4.50) 
gdzie 
- = aye + 6(cr-1)V + 2-5 گرگ بر‎ (4.51) 
fo = cr(cr«c)V9/ r - 6erv r I T * cr(3-2cr-c)y V T (44.52) 


Dwuwymiarowa temperatura krytyczna To Zwiyzana jest 


2 V następującym równaniem 


C 
138.3) 
-W/kT, = 1n(1 - —==——) (4.55) 
4¥ ,)1-۷ 22 | 
gdzie 
4 
اریہ‎ ee NON Pans © (4.54) 
K 2621-1 0 | 
q E ۱ 
-1)جه‎ V 7) 


Analiza równania (4.20) pod kytem ustalenia zależności 
liczby parametrów krytycznych od parametru niejedno- 


rodności r , prowadzi do następujących wyników 


(podanych w postaci tabeli) 
x Tabela T 
(0.2,1). (po “(2,c) 


przedział wartości cr (0,0h2') 


liczba 
parametrów krytycznych 


To 19 LIT 


W celu ilustracji, przedstawiono na rysunku 4.1 


, k. š 


ES ama 


e ter 
dan DIA 4 ۱ 


تسا 
-.—l‏ 


zo p 


wykresy funkcji f, oraz -f> dla czterech różnych 
parametrów r nalez&cych do wyróżnionych w təbi 
bı i dla c=b , Parametry krytyczne V. A 
-cw/kT. obliczone numerycznie z równania (4.50) 


dla różnych wartości (c) i (r) se następujące: 


Gab 


1 ۱ 
C 
0.99469 0.99126 26.78980 
0.94542 0.91396 13 .65325 
0.88724 0.83308 10.20121 
0.82692 0.76048 8.45277 
0.80656 0.73851 8.01300 
0.76611 0.69812 2.535196 
0.70713 0.64632 ` 6.58160 
0.652541 0.60408 6.04046 
0.61947 0.58045 5.75772 
0.60415 0.56980 5.65412 
0.58258 0.55516 5.46721 
0.52845 0.51904 5.06891 
0.50000 0.50000 4.86558 


Tabela II 


Liczba wartości V. wyznaczona numerycznie pokrywa się 
z liezbs podana w tabeli I, jednakże uwzględnienie 


warunku : 
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powcduje, ze otrzymuje się jedng wartość US dla 
danego parametru nieğednorodno$ci r (patrz tabela lI). “ 
Warunek (4.55) wynika z postaci funkcji F(V) - rów- 
nania (4.5). 

Na podstawie wyników Hoodia tabeli II 
można powiedzieć, ze wraz ze wzrostem heterosenicznoßci 
powierzchni (zmniejszenie wartości parametru r ) 
następuje obniżenie dwuwymiarowej temperatury kryty- 
cznej oraz podwyższenie wartości V. (limes Ves "BY 


r—) 
Rozważmy obecnie zachowanie się wielkości termo- 


dynamicznych w rejonie punktu krytycznego. Posłużymy sie 
4 w tym celu klasyczną teorię zjawisk krytycznych dla dwu- 
1 wymiarowych systemów, sformułowany po Taz pierwszy 

W 1973. roku przez u N 

Klasyczna teoria zjawisk krytycznych opiera się na 
rozwinięciu energii swobodnej helmnoltza w szereg 

i potezowy temperatury 1 ctęzenia powierzchniowego | ° 

Gest Ir. 1/۵ , 5 - pole vowierzchni adsorbentu ) 

Równania stanu korespondujyce 4 bragsa-Wiiliamsa و‎ 
Quasi-chemiczny i Hilla - de Boera aproksymac Ja 
zakładają, że ciśnienie powierzchniowe jest analityczną 
Tunkc ję neas I w 0 e krytycznym, zatem roziiwe 7 ۱ 
dest w tych przypadkach przedstawienie energii swobodne ل‎ 


Helmholtza w otoczeniu punktu krytycznego szeregiem 


Taylora: ł 
(4,56) ۱ 


i qe R 
og Ee + A à yi + zt 23 10 r Aq 9 T ) + 


" SH m A 8 A 
۵ q^ + -h, A Ke یت + او‎ E SE X Pi T el A 
2 一 , à 3 


za * p = 


Vs s 
۷ 
b 


= SDA 


W powyższym równaniu wprowadzono oznaczenia: 


À = A/S, D e fF- İŞ: کا‎ ۲ - 1, (4.57) 


oraz 


/ 0+07 ۱ : 


4: حب‎ ma (4.58) 
1) ۱ 1 j 
| or.” ör 
| 
| ( c - oznacza wartość w punkcie krytycznym) . 
W rozwinięciu (4.56) uwzględniono fakt , że 
Any = 450 = 0 (tj. warunek 4,47). 
Zajmijmy sie obliczeniem pojemności cieplnej Cp 
| dla krytycznej "izochory" | = dE 
Dla Tat, | 
23 ۱ ih e 5 
i A =z ^۵00 + A 4 + 75509 T + tes = Ag y (4.99) 
Stud 


. 2 نہ‎ d 
.. 1 9^ 4 ke Pul | (4.60) 
uu ` hi c 02 | 


Or 


II 


C 


Dla ot. w punkcie | = ۳ mamy dwie Lazy. 
W Lem przypadku energia swobodna Helmhoitza wyraża sig 


nDastepujsco: (^.61) ` 
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۱ I 
4 در‎ S | 
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2 y T ٦ ZU 
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eg 
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7 1 
dug: ( ES Te Jer! ) + o(| STI) (4,62) 


otrzymuje się 
2 
QA 9421 , 
y" ev eR A (4.65) 


Odległość między Ce 1 Or w punkcie nieciysiosci 


wyraża się zatem następująco: 


=. pxe 6421 | 
C. pa vh = IT è (4.64) i 


-. B2 - 


IV.5 Obszar temperatur subkrytycznyeh 3 ; 


2830 2007 , ze dwuwymiarowy GƏZ siatkowy może 
istnieé w dwóch fazach: rożcielczone j Jd Bu 
wanej. Dla danej temperatury T  ponizej temperatury 
„krytycznej Lo obie fazy 008 być w równowadze. Ciśnie- 
nie p. pray którym istnieje równowaga (ciśnienie 
gazowego adsorbatu) może być wyznaczone z izotermy 
adsorpcji przy użyciu dwuwymiarowegzo analogu reguły 
Maxwella równych port 2475 

MP 
in(pe) = | Tiet yar / (xg) 
Va i 
Edzie Ve dest wartością względnego pokrycia powierzchni 
V dla pep, po stronie "gazu" dwufazowego rejonu (po . n: 


 Czytek kondensacji) zaś V4, jest wartością V dla 


danego Dà po stronie "cieczy" (faza skondensowana) 
rejonu dwufazowego (koniec kondensacji) . 


Oczywiście 


BRO) = INP) = IR وص ا‎ = C a * (ase رہ‎ e 


Ostatnie równanie jest jednocześnie równaniem izotermy 
1 adsorpcji w obszarze kondensacji Vas V < Vr , która Zu il 
A 43 | ° e % > 
WM ten sposób rozpada się na tray części. Dla obszarów: 


| ۱ "2 . gy 
BB CV. oraz 1٦ L izotorma ads وک‎ 6 jz jest Gane 
A 


równaniem złusznym Cla obszaru temperatur c DEE TONA 0 
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Znajomość ۷۴ oraz V umozliwia znalezienie 
wyrażeń opisujących wielkości termodynamiczne 
3 w obszarze temperatur subkrytycznych. 

A W tym celu przedstavimy znane wyrazenia dla funkcji 


termodynamicznych w postaci dogodnej do zastosowania 


w obszarze kondensacji dwuwymiarowej. I tak, z relacji: 


i 1 | À s = i In ۱ `a 
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(tj. równanie prostej łączącej punkty AER i AN: 
na wykresie A(V) - V ). 
Wykonując różniczkowanie równania (4.65) wedXug wzoru 
o in(p,) 
Q = ×) کے‎ | (4.72) 
A Xa LM 


otrzymuje się wyrażenie na ciepło izosteryczne 


w obszarze kondensacji: 


M? 
| Qs اي‎ 
MMC | 
Qt = — s constans . . (4.73) 


Wielkość Ast można więc znaleźć przy użyciu metody 
` równych pól w ten sam sposób jak In(p,) : 
= Korzystając z równania (4.68) uzyskuje sie wzór 


š na entropig w następującej postaci: 
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| W obszarze kondensacği 
a (4.77) 


E, = À. + TS, = E(V,) - RON Me At. ere و‎ 


Wzór na pojemnośś cieplnu warstwy zaedsorbowane j 


otrzymuie sie z równania (4.76). Mianowicie 


(4.78) 
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Korzystając z równania (4,78) można wyrażenie (4.80) 


Zapisać w nieco inne i postaci: 
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Focstawiajyc powyższe równanie do równania (4.79). 
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Get Ne) = C (Vo) - Molstad) (4.85) 
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Uwzględniając ostatnie awa równania we wzorze (4.82) 
uzyskuje się prosty postać równania pojemności cieplnej 
w obszarze kondensacji ( równanie prostej przechodzącej 


| przez punkty Got) i CAVE) © 
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(4.85) 
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Otrzymujemy następujące wzory na Ve 1 — : 
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Vako końcowy wynik otrzymujemy: 
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ə BD = 


na wartości funkcji Cy(V) w punktach brzegowych 
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من‎ = Voy Vr, 


Korzystając z wyprowadzonego równania izotermy globalnej 
korespondujycego z exponencjalnym rozkładem energii 
na powierzchni heterogenicznej typu "random". można 


układ równań (4.65 , 4.66) zapisać wpostaci 
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uzyto nastepnie do obliczania wartogci funkcji 
termodynamicznych w obszarze kondensac ji dwuwymiarowé j. 
Wyniki obliczeü przedstawiono na wykresach ( rysunki 


4.3 - 4.6 ). 


IV.6 Wyniki obliczeń numerycznych . 


Rysunek 4.2.8. przedstawia typowe izotermy adsorpcji 
w temperaturach superkrytycznych. Można zauważyć 
znaczny wpływ heterogeniczności powierzchni w obszarze 
niskich stężeń adsorbatu przeğaviağqey się w szybkim 
wzroście adsorpcji wraz ze zmniejszeniem parametru 
niejednorodności r. 

Na rysunku 4.2.b. pokazany jest przebieg izosterycz- 
nych ciepeł adsorpcji . Widoczne jest, ze w miarę 
przesuwania się adsorpcji W stronę wyższych wartości 
wazzlednego pokrycia powierzchni znaczenie hetero- 
geniczności maleje, Zwiększa się zas rola oddziaływań 
między zaadsorbowanymi molekułami. Y rejonie pokryć 
powierzchni bliskich jedności oddziaływania te 
całkowicie decyduj O przebiegu zjawiska. 


Rysunek 4.2.c. ilustruje “Pty heterogenieznosci 


— powierzchni na pojemność cieplny warstwy zaadsorbowene y. 


zlędniające wkładu pochodzącezo 
oh trójwymiarowych oscylacji zaadsorbowanygon 100 1664 . 
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wynikający z równoczesnego istnienia oddziaływań 
między zaadsorbowanymi molekułami i heterogeniczności 
powierzchni. Pojemności te nie —QÓ również T. 
wkładu pochodzącego z uvzglednienia perturbacji 

pionowych wibracji molekuł, spowodowanych różną 

wielkością minimów potencjału adsorpcyjnego na po % 
wierzchni adsorbentu. Ten typ wkładu obliczany był 

już w literaturze przez Draina i Morrisona. 

Typ krzywych przedstawionych na rysunku ۰ 

tj. pojemności cieplne Gud Mk wykazujące wyraźne | 1 
maximum , 2 obserwowany w przypadku adsorpcji He” 1 
na wstępnie pokrytej argonem niedzi/ ^0 , 

Rysunek 4.2 ilustruje wpływ siły wzajemnych SSES ۴ 


ływań w warstwie adsorpcyjnej na wielkość adsorpcji (a), 


ciepło izosteryczne (b) i pojemność cieplną (c). 


Wyraźnie widać przewagę sił wzajemnego oddziaływania 


w rejonie średnich 1 wyższych pokryć powierzchni 4 
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nad efektem heterogeniczności. W momencie gdy niejedno- . | ۱ 
rodno$6 powierzchni traci znaczenie „wykresy funkcji. 

| REN termodynamicznych upodabniają sie kształ tem. d 
do wykresów dla powierzchni homogenicänych. 
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tycznych pokazany jest na rysunku 4.4. Interesujące 
_ jest, że ze wzrostem heterogeniczności powierzchni . 
— obserwujemy szybszy liniowy wzrost pojemności 

- cieplnej w zależności od pokrycia powierzchni 

(w rozważanym obszarze kondensacji). ] 
Na rysunku 4.5 pokazano diagram fazowy gaz - ciecz 

dla powierzchni heterozenicznej scharakteryzowane j 

parametrem niejednorodności 200 . 

Na następnym schematycznym rysunku 4,6 pokazano 

jak krzywa koegzystencji zmienia swój kształt 

T wraz ze zmianę stopnia heterogeniczno$ci powierzchni. 

© Ostatnie wyniki (rys-4.5, 4.6) wskazują na to, 

że heterogeniczność powierzchni niewiele zmienia 
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V. TROJWYMIAROWY MODEL ADSORPCJI. 
ADSORPCJA MOBILNA NA  HETZROGENICZNYCH 
POWIERZCHNIACH. 


V.1 Sformutowanie tro jwymiarowego modelu adsorpe ji 


Dotychczas teoria adsorpcji na powierzchniach 
heterogenicznych bazowała na dwuwymiarowym modelu 
adsorpcji, gdzie izoterma adsorpcji była wyznaczana 
z warunku równowagi między faza objętościową a fazą 
 zaadsorbowaną à 
TW trójwymiarowym modelu adsorpcji, przyjmujemy, 
Ze cząsteczki gazu mogu poruszać się w trzech wymiarach.- 
Tym samym, opis takiego układu staje się trudniejszy 

niż opis modelu dwuwymiorowego, opiera się on także 
na zupełnie innym formaliźmie. Zasadniczą charakterystykę 
Układu adsorpcy inego W trójwymiarowym modelu adsorpcji 


tworzy jedno i dwucząsteczkowe funkcje rozkładu n (y), | 


1 


8620). ) . Znajomość tych funkcji pozwala na wyznaczenie — 
1 u | 


“M‏ ے 


Wszystkich charakterystyk termodynamicznych układu, 
tak; izoterma udsorpcji wyraża się przy pomocy 


SE ocząsteczkowćj funkcji rozkładu (gęstości lokalnej) 


m dée, 
Tównaniem/ V. | NH. TN ot : 
"SAB ۱ eat جار‎ - n)dr ۱ (5.1) 
r / سے‎ E 
L V " 

bl ۱ 1 

Bie N oznacza ilość zaadsorhowaną, | n. "nn Fe 


E » ۱ 


PM m. e ° ۰ = LÀ s 5 4 = 7 ون‎ 7 u > m ia r 3 
BeStosciy fazy objętościowej 263 E Jest rek torem E 
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położenia, 
średnia energia wewnętrzna” 1/ 
(5.2) — 


E = 2 RT wi n (p)v(r)ar + RZE Eo) u(r45)àr4ür, 


u ۳‏ — .== > سے 
b -‏ 


gdzie v(r) - potencjał działający na molekuły układu, 
u(r; 3) - potencjał EE onte dwóch وا‎ 

molekuł. 5 = tę Z Dal 

Sprecyzujmy dokładnie ER trójwymiarowy. Zakładamy, że/ ^4 

a) ciało stałe tworzy Sztywną, nieprzenikliwą dla 
cząsteczek gazu strukturę, oddziaływując na 
cząsteczki gazu potencjałem v(r) , który musi 
mieć tym samym charakter potencjału ze SC STE 
ścianą (interesuje nas 3 n ls. aà nie absorpcja), 

b) cząsteczki gazu sy kuliście symetryczne, a więc 
energia oddziaływania par uri) zalezy tylko - | 
od odległości cząstek, | bag ۱ ` 

c) energia potencjalna N cząstek moze być przedsta- 


A 


wiona w postaci Sumy 


i 

>] R 3 ER 

g MA 

A APT 

działające sity S4 si Xani. FOR EVE 

d) ws szystkie Ch ۱ 
m namiczne à = 

zasiegu, a wiec wielkości termody: 2 x 

ıb E 

nie zależy od ozraniczeń przestrzennych układu. ` K 
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Jeśli 1 zazu ss bardzo niewiel! temper ol 


Beary dostatecznie wysokie TO gęstość lokalna moze by i 2 
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przedstawiona równaniem Roltzmana/ 1/ 
(1) | | | ; 
n` ریا‎ = .n-exp(-v(r)/KT) (5.4) 


Pray wyższych gęstościach gazu odchylenie od prawa 
Boltzmana staje się znaczne, w tym przypadku gęstość 


| lokalny można zapisać w postaci szeregu wirialnego/ ^» M 


nV.) P g(r, (ary 1D Co, رھوں‎ qe) 
m=2 
gdzie 
| "exp(-v(r,)/KT) Tr; ۷ | 
mec) = |. r gi . (3.6) 


V jest objetoécis gazu zamkniętego przez płaszczyzny 
ciała stałego. 
Praktycznie, stosowalność równania (5.5) jest ograniczona 


do obszaru niewielkich gęstości gazu. Wyznaczenie gęstości 


lokalnej dla dużych gęstości gazu i w niskich tempera- 


turach jest skomplikowane i może być jedynie uzyskana 


T "R m 


w prostych przypadkach homogenicznych adsorbentów. 


Z tezo względu, zastosowanie równania (5.1) jest 


| Ozraniczone do przypadku powierzchni homogenicznych. 
| W rozdziale V.2 przedstawione zostaną krótko sposoby 
| 


| wyznaczania rozkładu gęstości lokalnej dla powierzchni 


| homogenicznej a następnie zostanie skonstruowana metoda 


tad ści powierzchni hetero- 
wyznaczania rozkiadu gęstości dla powierzchni netero 


genicznych i w konsekwencji izoterm susorpe jt 


1 — " 7 "vy ( v ' روہ‎ vo? ) 1 
na powierzchniach nsterogenicznych | POL ODL 27 
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A Pokazane zostanie również. przejście od równań 
A, uzyskanych dla modelu trójwymiarowego do całkowego 


równania izotermy adsorpcji. 


| V.2 Gęstość lokalna dla powierzchni homogenicznych 


Gęstość lokalna dla powierzchni homogenicznych 
D moze byé wyznaczona przy uzyciu następujących teorii: 


a) Born-Green-Yvon (BGY) / 5-8/ 


b) Percus-Yevick (PY) 

c) Hypernetted Chain (HNC) /9-12/ 
W kazdeğ z tych teorii zagadnienie wyznaczania gęstości. 
lokalnej sprowadza się do rozwiązania b... 


numerycznego) odpowiednich röwnan całkowych. I tak, pa y/ Di 


Die Je (e Mr, T5 ini? (x4) J Cu 212)/«D dry 9‏ رررع) P‏ مھ( 


V رخ 7) ب‎ )/k?) (5.7) 
gdzie gs” (z. 125) - dwuczysteczkowa funkcja korelacji AE © 


V4 - operator rözniczkovania po wektorze lu” 


/ * gezu 

potozenia 24 ٠ 

ü f gi ğ 

W celu rozwiązania powyższego gna zastępuje się Mo E 


AUE: E 


D 
GIE przez radialna funkcję dystrybucji gt? ie 
ein 


objętościowego płynu. W tej aproksymac ji 8 بی‎ % 2 


Jest następujące: ; 
( ) 

کر کی E‏ | 

a V4ln(n (r4) - İst” anon ) Cr Jexp( لل‎ E 


a 


m | 1 
> oNQG(GAM KP) ٦ ee کر‎ 
- E R 7 Ê | 


— 


ə ہہت‎ 


mür kü a ۳:23 Wire 1x ور‎ ! 


bat‏ | ہریڈ ga‏ ری سر Zeg,‏ ود 


d 4 | ind Z Wa Su nö er | "Pet, 
LT 3 We E 2. gu 


d | BOJE 


1 Przyblizenie użyte powyżej polega., - ozólnie rzecz 
biorąc, na zaniedbaniu wpływu 0 St iney gaz- ciało, 
4 state na 'korelac je dwuczasteczkowe. Wy Jasnijmy to na 
przykładzie pierwszego współczynnika wirialnego in p © 
w rozwinieciu dvuczasteczkoveğ funkcji korelacji. ظ‎ 
W tym przypadku, współczynnik ten może być zapisany + 


następującą de rə 


| (Ta )f(r44)f(r54)02, وی‎ 
gdzie f(x; 5) oznaczają funkcje Mayera DEL 
d | EN 
f(r,.) exp(-u(r,.)/KkT) - 1 (5. 10 E 
=1 j lJ PER 
۰ L. L* gt 
5 pam ١ 


Zastapienie dwucząsteczkowej funkcji korelacji grana) ۰ 


1 
| 


e 
przez radialny funkcję dystrybucji ge (245) 


jest równoważne zastąpieniu całki (8. 2) przez nastqpu- 


jęcą całkę 站 wat: RE 
a 
f(r4s)f(r» ور‎ (12 : x 
3 ۱ ba? 
š Jak widać, różnica polega na braku w rów. (end MUN 


na nieuwz;sledn). eniu oddziaływań gaz- WETE — E 


a 


WE? 


| (czyli 
; nia. W calce ostatniej 


“M oraz na zmianie granic caikowa 


granice te ae ty 0 od zasięgu oddziaływań nędzy” 


| czę ‚steczkowych (symbol 00 x natomiast — 


A KA bd 
wchodzę nam nowierzchnie ciału stałego. EE 


Powyższe sproksymac ja jest więC یئ‎ bardziej słuszna 
"s - 
m. 
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im słabsze są oddziaływania gaz-ciało stałe (niższa 
energia udsorpcji lub wyższe temperatury). 


Teorie PY i HNC bazują na rozwiązaniu 


następującego układu równań/ TV, ۱ 
rownania Ornsteina 
. | some) 
h(z) = c(z) + " h(z4)c^" "Lëns dx. (5.12) 
oraz 
h(z) = c(z) + (h(z)«1)exp(v(z)/kT) - 1 (PY) (5.13) 
h(z) = e(z) + la((h(z)+1)exp(v(z)/kT)) (HNC) (5.14) 
gdzie h(z) = PISTE i 
z - odległość od powierzchni adsorbentu 
c (pau) - pez poé retina) Turkce JAM Kkore TE 
dln gazu objętościowego (znana z teorii faz SW 
ob 0ڑ‎ oácioxych). r, wr 
Moria pr^ FNE przy judy nastepujycy model. ` m d 
| układu adsorpcyjnego: Układ edsorpey jay desə” pw Paa 


L ۰ ۰ de ۰ ] " P al y | m à 
rozpatrywany jaks “şuskxadnikova mieszanina, przy czym. 
İ 5 r * 4 


jeden skXadnik jest reprezenzowany przez pojedyfczą 


d nd $ pe - WW ETE rb di mm 
5 dużą caystke. Rozważając ET 4 cay promień tej JE 


astî döy do nieskończoności uzyskuje się układ 


L< powierzchnią sraniczny. ^ ten وی نوک‎ p اف چا دہ‎ 
3 zt uet m? s ES 2 
i UNG ala homogenicanych miesaunim, użyte PEL a 


Tw grenieznym przypadku wyżej zdefiniowanym, dają nam | 
Możliwość opisu تل‎ 77 udsorpey ان‎ CR. eebe ۳ 
N ee "pil AC ١ ا ہے‎ 
EUM uS, nadniezić, 2€ rozkrad gestofci loxsineg uzyskany + 


" 1 ° " WT" - 
a é x 1 | > KE 4 i ~ a . 
E ا ره‎ f, Li A " S am ^7 zə à 4 ) ! 
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jest w najlepszej zgodnośći z wynikami uzyskanymi 

P na drodze symulacji kemputerowych. Oczywiście, jak 
zostało wspomniane wcześniej, ścisłe rozwiązania 
układów równań PY oraz HNC można otrzymać dale 
w przypadku homogenicznych adsorbentów. Czytelników 
zainteresówanych metodami numerycznymi rozwiązywania 


tych równań odsyłamy do prac orginalnych 919, 


V.3 Gestoás lokalna dla powierzchni heterogeniczne j 


Załóżmy, że całkowity potencjał adsorpcyjny v(r) 
/ 14/ 7 


można przedstawić w postaci 


v(r) = v (zə * v4Cr) (5.15) 


gdzie vo C2) jest odpowiednio wybraną funkcją 
dla homogeniczne(50 układu odniesienia, 


və) jest tzw. heterogenicznym Wd BERN 


1 Takie traktowanie oddziaływań gaz-ciuio dares jest 
analogiczne do tego, które wystepuje w teoril gazów 
و بس ہت‎ NS Główna różnica między " perturbacy jnym" 
podejściem zastosovanym tutaj a prezentowanym 
x 0 w pracy /15/ 8 na tym, że w układach gazowych 
podział całkowitego potencjału na część "odnoszącą" 
i zaburzającą jest 6 określony, natomiast dla 
A układów gaz-ciazo stałe podział taki moze byś użyty 
| Jedynie dla powierzchni krystalicznych. 
Je Dla powierzchni z nieperiodyczaym potencjałen istnie je 
więc konieczność skonstruowania nowej metody podziału 


potencjału lub z;módyfikowenia netody sveele >a 
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Pepe, Fa پو‎ 5 nə. Kai ` e n 
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Jeśli stosujemy zgodnie z równaniem (5.15) potencjał 
zaburzający va4Cr) w homogenicznym układzie 

z potencjałem رم‎ wówczas całkowity potencjał V(r); 
zmienia się od v C2) do vg C22 ev4Cr) iw Kan rri gr 
gestoóó lokalna zmienia siq od nC 7 ری‎ do nC Ur). 


"bar 1 | 
Rozwi jając m. ( (1) “expiv, (2/0) w szereg funkcjonalny 


Taylora względem exp(-v4(r)/kT)-1 i używając równań 


AC 
ol wy prowadzonych przez Percusa/ 1% | otrzynujeny 


alte) expr (/ = ng (2) + (5.16) 


1 | 
1 141 / 

> TT p 1 [ Cexp(-v4 Ce ,)/ KT) - 1)dr4...dr; 

i=1 jz=1 | ۱ 


gdzie p sa tZW. i-czasteczkowyni funkcjami 


rozdziału Urde11la 1617. 


Dla naszego ukiadu 


p (n4. oz) A | €V 
R دوه وال‎ 
p. j (5.17) 
O daag ۳0 re | 


J 


Równanie (5.16) może służyć do znalezienia gęstości 
lokalnej ۰ج‎ 12 (r) dla powierzchni heterogeniczne]j. 
W przypadku powierzchni homogenicznej dostaje się: 

n (y) = یس‎ zzolnie Z przewidywaniami. 

Człon równania (5.16) zawierający i-czasteczkowg 
funkcję Ursella pg, reprezentuje korelacje między 


wO = 


(i) cząstkami i wpiyw niejednorodności. powierzchni ` 


telas r 33 £2. 


v 


| | = e = ? 


na te korelacje. 


Dla wystarczająco małych wartości potencjału v(x) 


(mała niejednorodność) oraz dość kysoki ch tempóratut, 
szeres (5.16) jest szybko zbieżny, możliwe jest zatem 
w tym przypadku zaniedbanie wszystkich członów szeregu 
z wyjątkiem pierwszego. Postępujyc w ten sposób 


otrzymujemy : 


2 


nC (y) = exp(-v4Cr)/ KT) PCE) (pp) (8xpl-v4(1,)/kT)-1)dr, 
n(n) est eating) (5.18) 


Powyższe równanie można jeszcze bardziej uprościć, 


zakładając, 28 dwuczysteczkowa funkcja rozdziału pt2) 


187, 


wyraża sią naàstqpujuco 


A e ` 7 e. $ | 1 a 
Ty Ed c c “(İe-r-l کربت‎ D Vo 


i ; 2ھ‎ 
AY جا‎ 

2 ۱ LE m A 
(2) jest radialna funkcją dystrybueji w fazie ای ۔‎ 


gdzte Ly 


objętościowej. x SEE 


Podstawiajnc równanie (5.18) do rómania (5.1) otrzymuje 


Sie równanie izotermy easetpegs * [ ۱ 


AD 


ANIE RIA ADO ją 3 


= 1)(exp(- te YE) ars) - 2 (5.20) ` 
| l F 


۱ e AA 1 Sal 3 MET r ا‎ ə 
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powierzchni ciała stałego. 


W tym przypadku równanie (5.18) przechodzi w zwyczajny 


rozkład Boltzmana 
(1) A E | 
n` ^(r) = n"exp(-v(z)/kT) (5:2 


W konsekwencji izoterma adsorpcji wyraża się 


następująco: 


N = a İ (exp(-v(2)/4D)- az (5.22) bs 
: | 
Przyjmując dodatkowo założenia, że 
(a) miejsce adsorpcy jae jest ou IR 


przez minimum energii potencjalnej wzgledem 


ruchu prostopadiego do powierzchni, 
tj. minimum v(2)szB رہ‎ 

(b) miejsca adsorpcyjne Uważa się za Heme samego 
rodzaju jesli okreslone su przez ty samą 
wartość E, 1 باه‎ 


mozna vównanie (5.22) zapisać w postaci: ` 


ə . ch NA. 
która jest jak widaś izoterma eo eh LE deer BUT 
heterogenicznych! | : x | 1 r At 

(B) Brak korelacji niędzycząsteczkowych. Sai | ix 

3 W tym przypadku ski 2 =0 i równanie (5.18) przyjmuje ۱ 2. 


۲ Pos tac 


TOR - و‎ ) (z) ergi )را‎ 0/1 
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Izoterma adsorpcji (5.20) upraszcza sie obecnie 


` do równania: 


È (nf co exor (an - mar (0 (5.25) 


Przyjmując założenia (a-b) z punktu (A) oraz 
zaniedbując adsorpcję vynikaduca z istnienia potencjału 
"sztywnej ściany", równanie (5.25) można sprowadzić 


do następującej postaci: 


N = | £(n,K, (E) E, EAE) (5.26), 


zi ds esst, وو وی‎ - Vax (5.29) 
V 
jest state Henryego , 


司 2 = 5 C d | 
nd( 1-b)K,/ CL Cot 1) 7 . 5 (528) 


4 


一 
p 
PN 

də 

مح و 
em?‏ 
Ko‏ 
RNA‏ 
Na‏ 

D 


natomiast d oraz 5 oznaczają ocpowiednio średnicę 
cząstki adsorbatu i pole powierzchai adsorbentu. 
` t 2 b. G 


1 ` > re 7 z. [ts ط2‎ 
Calkovu reprezentecoç izotermy Şiobalned (5.36) 


można uwzzaé za analog równania er Aach ho 7 
Stud tez mozemy powiadzieć, Be réwnanle izotermy ۱ 
&lobalnej (2.1) moze być uży ta JE کا کو‎ opisu | 
Uuktadéw adserpcyjnych w których əfokt adsorpe ji. 
Wynikagjacy 2 igtnienis odd zia; wena prayolugs ines > © 

| powierzchni jest dużo Ya ود عن‎ nig afekt, وس‎ mj fia 
zwiazany Z potent jalen 8% 2 +11 6 ہام رہ‎ onen done ۹ 


ZAŁ > 
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- SZ 
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k 
4 


EZ s 


mieduy gazowymi molekuzami są zaniedbywalne, 

Poniewaz ostatni warunek jest spełniony wówczas, gdy 
oddziaływania między zasdsorbowanymi. molekułami sy mate, 
posługiwanie się równaniem calkowym (2.1) lub (5.26) 
do badania wpływu oddziaływań udu A Ean 

na całkowitą adsorpcję, moze prowadzió do przyblizonych 
wyników. 

Równanie separacji potencjałów (5.15) nasuwa również 

myśl traktowania całej części yprzyciasadacel potencjału 
jako perturbacji, 5 układem niezaburzonym byłby 

układ ze sztywną 01۰ Zagadnienie to, było w sposób 
fragmentaryczny rozważane przez Kuniego i Rusanola/, SE Wz 


Uzyskali oni następujące wyrażenie na gęstość lokalną: 
| P, 
n(1)(2) = n - n^ v4(2) (5.29) 


gdzie e jest izotermicznym współczynnikiem 4 


AL.» 


x =. 12+ 一 
Y) JP T - 
: 3 E / * ^4 نم‎ 3 2 342913 e : DG f. 
Równanie (5,29) geet podobne do wyrażenia Qaxıe otrzymar rf 
: 4 inzanie równanie na-Greena- £ ch $ 
ا‎ * przez rozwiązanie równania. Borna-Grgena Yvona, MK 


ia vestoóci lokalnych i izoterm | 
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v(r) = v,(2) + v4(r) (5.15) 
OO dla 2¢0 i 
v (3) = " 2 3 : (59750219 
t zz z» 0 
A 
var) = E ZA (a eos (ERES) * cos( STZ) (5.30.2) 
K y ` 


Fotenc jaz ۷ (2) jest potencjałem Sutherlända, natomiast 


postać zaburzenia və (z). jest podobna do potencjału 


Ą 23 Y š 了 1 ` : 

uzyvanezo przez Hilla’ / Go modelowania powierzchni RE ir 
yd A zy 
z periodyczny zmianq pov tenc ۰ Stałe 5 i y iz | i 
مم‎ 200 = 

sieciowych) , Ë jest parametrem zwinzanym 2 WE cd 


potencjału zaburzajycego. Maksymalne wahania aem =, 
v(r) wynoszą więc A (a,+1). í Czynnik Gy (+2, y” A 
powoduje szybki 8 ۷ potencjatu zaburzającego Ei 


ze wzrostem odległości. od powierzchni - szybszy niz  ăč 


d € Ae ON ` š P ES 
Spadek po vente jału v C2) e Parametr € W اس نیک‎ 


potencjału v (z) jest ۵ energiy adsorpcji 
5 0 W 


na powierzchni 2 potencjałem vir). 


Jako potencjat oddaziazywania u(r) dwóch gazowych = 


molekuł przy 0 potenc jet sztywnych kul. 
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genicznej scharakteryzowane j potenc jałem ۷ (2) oraz 
radialnej funkcji dystrybucji g(r) . Ograniczając 
badania do obszaru bardzo małych gęstości adsorbatu 
(1077-107?cząstek/A”), możemy radialną funkcję dystry- 
bucji dla sztywnych kul przedstawić w postaci 3 


eC?) (y) = 1 + g$2 (r)n + ES (rn? (5.50.4) 
gdzie gi”) są powszechnie znanymi współczynnikami 
| بت‎ R 


2 
Jakkolwiek wartość EN er) 


w zakresie rozważanych 

gęstości jest bardzo bliska jedności, nie możemy przyjąć 

tej wartości jako przybliżenia ponieważ we wzorze (5.18) 
mnożymy ją przez bardzo dużą wartość funkcji exp(-v4/kT). 
Gęstość lokalną dla powierzchni homogenicznej wyznaczono 
przez numeryczne rozwiązanie równań PY (5.12, SMS”) 


metoda iteracji stosujac jako zerowe przyblizenie 


pierwsze cztery człony rozwinięcia wirialnego: 


v2) 


n 2) - exp(--2 )(n + 2 (z)n* * 345 ۱ (2)n? + 


< ” 


: 2. ın. (5.30.5) 


zdzie aC 2) SĄ współczynnikami virialnymi,danymi 
i 


za pomocą następujących całek: 


21w51)(z) = |)6)22(-1(۶)22( 422 | 5 


— 


- 


e 


DKG (g(z4)-10| glzo)f(112)F(Ta3)dT2dTz + 
CA ə 
+. mE | oa? -1 (1 (2,5) (£03) f(z44) 8750 rs 

00 | s سس‎ E ۱ ۱ MR. dnd 


< 4 r a e I b. ra 4 
EA z 1 نے‎ A سی"‎ Y = Fs mo. 2 
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Zapisanie wszystkich całek wchodzących do WYD (z) 
zajęłoby co najmniej dwie strony, można je znaleźć 
۷ literaturze/ 2: 12M | “R” 


W powyzszych rövnaniach 


| (-4 dla ef Ko 
| fir, ;) = ijt 
Ë d Tij) 1 


1 1 ët: 
Gałki we 2 us 7 "c ) obliczono numerycznie stosujac 
metode iteracji Gaussa. | 


Przez całkowanie gęstości lokalnych dla powierzchni 


homogenicznej i heterogenicznej otrzymujemy zzodnie 
z równaniem (5.1) odpowiadające im izotermy adsorpcji, 
W obliczeniach stosowano następujące stałe para- 
| metry: 8421.0, © 0 cal/mol, z =1.0, T=100 K, 
a (wartość parametröv:(z.) i (a) podana jest 
w ğednostach średnic molekuł udsorbatu ). Jako przedziar 
zmienności x i y , wystarczy wzięść przedział (0,a/2) 
co wynika z okresowości i symetrii funkcji CAE) 
względem zmiennych x i y. | 
b Przyjmujgc, ze parametr £ =3000 cal/mol, wyznaczono 
rozkłady gęstości lokalnej dla powierzchni homogenicznej 
oraz heterogenicznej i porównano je z rozkładem gęstości 
gazu doskonałego (rozkład Boltzmana). 


I tak, rysunek 5.1.4. przedstawia rozkłady gęstości 


lokalnej dla gęstości adsorbatu n-407 8 یو‎ tek/A? 


(odpowiada to ciğnieniu sl Tr w T=100 K), wzdłuż . 
wyróżnionych prostych na powierzchni adsorbentu 


(x/a=s0,1, 0.3, 0,5) W bezpodredaim sąsiedztwie powierz- 


= 111 > 


chni (220.0). . 

Rysunek 5.1b ae" analogiczne krzywe rozkładu , 
gęstości lokalnej dla innej wartości zęstości | 
adsorbatu- n = 107 Tcząstek/A”. 

Rysunek 5.1c pokazuje zależność rozkładów gęstości 
lokalnej od odległości od powierzchni adsorbentu. 
Obliczenia wykonano m——- następujące wartości 


Ocząstek/A”, x/ 920.055. 290.102.2055 


parametrów: 5-10 
Rysunek 5.2 przedstawia okresowe zmiany gestoéci 
lokelnych obliczonych 4 0-1032 Pezastek/À^, x/a=0.3, 
z-0.0 i trzech różnych wartości stałej sieciowej a: 

a= 1, 2, 10. Frzebieg krzywych pokazany na tym rysunku 
jest odbiciem okresowych zmian potencjału v(r). 

Zgodnie z przewidywaniami, 29 wzrostem ęstości ad3orbetu 
(wzrostem ciśnienia) zeştoğci lokaine wykazują rosnące 


odstępstwa od rozkładu Boltzmana. Można takže zuuwaźżyć, 


że odstępstwa te sy większe dla większych wartości 
potencjału GK Q Fakt ten tłumaczymy wzrostem 

2: ۱ : 

oddziaływać miedzycz&steczkowych na skutek zagęszczenia 

E. molekur w obszarach dupy cb wartości potencjału VCE). 


W miare wzrostu odległości od powierzchni adsorbentu, 


zanika wpływ zmien „otencjału v(r) i tym samym 


rozkłady zęstości lokalnej upodabniają siq swoim 
kształtem do rozkładu gęstości 0 doskonałego. 
Rysunek 5.10 pokazuje jak bardzo wartości gestosci 
lokalnej zalez% od odlegtoáci oc powierzchni. Zmiena 
wartości parumetru (2). od 0.0. do. 0.2 + powouuje 


zmnie szenie wartości gostoáci lozalnej około 10000 
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razy. Zależność uoa z pokazana na tym rysunku 
jest charakterystyczna dla adsorpcji monowarstwowej - 
jedno maximum w 22000 . ۱ l 
1 Obliczenia izoterm adsorpcji przeprowadzono używając 
różnych wartości parametru £ (.i.=2000, 2500, 3000 
۲ mal). Dla każdej wartości ۶ obliczono: izoterme 
Henryego (koresponduje ona z rozkładem Boltzmana), 
izotermę dla powierzchni homogenicznej o energii 
adgorpcji równej £ i izotermy dla powierzchni hetero- 
| genicznych różniących się wartością parametru a. 
i Przebieg tych izoterm pokazany jest na rysunkach 
5.5a,b,c. Zgodnie z oczekiwaniem, adsorpcja na "povierz-. 
n chniach heterogenicznych jest wieksza niz na powierzchni 
homogenicznej o tej samej średniej energii adsorpcji. 
| Zwiększenie okresu potencjału (stałej sieciowej a) 

| powoduje zmniejszenie adsorpcji na powierzchni hetero- 

| genicznej. Fakt ten można by tłumaczyć w modelu dvu-. 
wymiarowym adsorpcji, zwiększeniem odległości między 
centrami adsorpcyjnymi (minima lokalne potencjału 

*  v(x,y,2) ) - zmniejszeniem oddziaływań bocznych 

w warstwie zaadsorbowanej. 

5 Wyjaśnienie tego w oparciu o 0-581 trójwymiarowy 
adsorpcji można przeprowadzić posługując się schematycz- 
| nymi rysunkami 5.^ Í 5.5 (nastepna strona). 

Rysunek 5.4 przedstawia (fragment rys. 5.1) „liczbę 
zaadsorbovanch molekuł znajdujący sie w zasięgu 
oddziaływań (na rys.zakropkowane koło) wyróżnionych 


molekuł (czarne punkty.) » Widoczne jest, że w przypadku 
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a=2 obszar oddziaływań obejmuje większą „liczbę cząstek 


zaadsorbowanych niż w przypadku gdy wartość a= 10. 


Tak więc stała sieciowa ma duży wpływna całkowitą 
energię oddziaływań bocznych a tym samym i na proces 


adsorpcji. 


Rysunek 5.5 pokazuje obszar energetyczny powierzchni ` 

r leżący w zasięgu oddziaływań zaadsorbowanej molekuły 

| (czarny punkt). Z rysunku tego wynika, ze dla dużej 
wartości a, obszar ten obejmuje jedynie miejsca 

|  adsorpcyjne o zbliżonej energii adsorpcji, natomiast 
dla małych wartości parametru (a) może on nawet obe j- 

| mować miejsca adsorpcyJne o wszysy£ich możliwych 


wartościach energii adsorpcji. Tak więc, przypadek 


pierwszy odpowiada wymogom płatowego modelu tópo- 
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graficznego powierzchni heterogeniczneğ (płaty jako 
zbiór miejsc adsorpcyjnych o zbliżonej energii adsorcji). 
| Przypadek drugi (duża wartość parametru a) odpowiada 236, 
| modelowi o pośredniej topografii miejse adsorpey jnyéh, 
którego skrajnym przypadkiem jest modeli"random" 
powierzchni heterogeniczne j. Zatem fakt wzrostu adsorpcji 
ze wzrostem wartości stałej sieciowej (a) może być 
| tłumaczony różną topografią powierzchni heterogenicznej. 
Rysunek 5.3 pokazuje również, że wraz ze zmniejszeniem 
średniej energii adsorpcji (podwyższeniem temperatury), 
różnice między izotermami adsorpcji obliczonymi dla 
różnych wartości parametru (a), zanikają. Jest to 
oczywiście spowodowane zmniejszeniem ilości zaadsorbo- 
wanej a w konsekwencji zmniejszeniem udziału oddziaływań 
w fazie zaadsorbowanej. Zatem w tych warunkach, parametr 
(a) a więc i topografia powierzchni heterogeniczne]j 


ma niewielkie znaczenie. 


UREA 


V.5 Drugi współczynnik wirialny dla dwucentrowej 
powierzchni heterogenicznej 


W rozdziale tym zajmiemy sie zagadnieniem 
wpływu oddziaływań cząstek zaadsorbowanych na różnego 
rodzaju centrach adsorpcyjnych powierzchni hetero- x 
genicznej na wartość adsorpcji fizycznej. 

Punktem startowym do naszych rozważań jest rozwinięcie 
wirialne dla izotermy adsorpcji, dane następującym 


równania 
2" iW Ax 


—— (5.31) 
J i 
xor. 
J= 


N/s = [”= Kun + 


gdzie N - Srednia liczba zaadsorbowanych ozübtək, 
S - pole povierzchni adsorbentu 
(7- stężenie powierzchniowe (adsorpcja Gibbsa) 
K,- stała Henry'go w trójwymiarowym modelu 
adsorpcji | 
B.- nieredukowalne całki Mayera 
n - gęstość fazy objętościowej 
W.- współczynniki wiriału gaz-ciało stałe. 
Zakładając, że energia oddziaływania gaz-ciało stałe 


dane jest równaniem 


00 z 4 0 ur 
x € (2id,2id+d) 


3 d 
| Dep I 0,h) (5.52) 
v(r) zei x € (2id+d,2id+2d) Í 
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gdzie 1-0, 11 ^ Pda, 


mamy zgodnie z równaniem (5.6) 


EN dla gx 0 


8 x E(2id,2id+d)  ” 

g(r) = (4| z €(0,h) (5.32) 
STI x £€(2id+d,2id+2d) 
1 zyh 


Równanie (5.32) określa jednocześnie topografię 


powierzchni heterogenicznej, której model jest 


| 


przedstawiony poniżej 


- 


Rys. 5.6. Model powierzchni heterozeniczne), określony 
przez potencjał oddziaływania gaz-ciatlo stałe 


(równanie 5.32) 


Powierzchnię taką można uważać Za zbiór homogenicznych 
płatów określonych przez głębokość studni potencjału 

(ves vyr): Gzerokość studni h jest mniejsza od jedności 
(w jednostkach średnic cząstki adsorbatu) jeśli zakłada 
się, że zasięg six powierzchniowych nie wykracza poza 
monowarstwę czystek zaadsorbovanych. 


Niech potenc jaz oddziaływania gaz-gaz będzie potencjałem 
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4 سی وس 
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. نے‎ 20 RA (e 
k əzan FER > 

۱ =" 4 
5 t , en, 


AT ME 


sztywnych kul 


oo dla ۲2 د‎ 1 
R 12 
ukz12) = {> : | (5.54) 
| Taka postać potencjału U(r45) określa nam waruńek d>2, 
oznaczający możliwość kontaktu molekuł zaadsorbowanych 
tylko na najbliżej susiadujucych centrach różnego 


i rodzaju. Zgodnie z równaniem (5.10) funkcia Mayera 


-1 dla 1۳12 < 1 


AE PEO í 
12 | 0 RE (5.35) 


Zağmiery się obecnie obliczeniəm całki Wo - drugim 
współczyanikiem yiriażu s na OTE EVE ttniemy-rożzwinisgoie 
wirialne izotermy adsorpcji (5.31) . Definiuje go" 


następujące równanie 


0 | | 
A | | 
gdzie 
OSZE 2775 | 0 
co 
€ Symbol CO pod znakiem Cułxi oznucz8, ze catkowanie "i 
d dest rozciagniete na cal: praesti/ef konfiguracyjny. 
Analiza równa potencjałów wykazuje, ze całka Ho 
jest Plana od zera w następującym obszarze całkowania : 
j X4 G Ax X>.€ (x.-1,x481) 
y > À y Yo& (Y4- 5334 +1) | (5.58) 
z, € (0,h+1) zə &(2,-1,h) 


Ponieważ wyrażenie podczakowe (rów. 5.36) jest funkcją 
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periodyczną zmiennej Zu O okresie równym ed و‎ 

| wystarczy rozważać granice całkowania dla Xa przedziar 
(0,24). Dla x4 € (0:4) można na podstawie równań (5.33) 
oraz (5.35) zapisać 


F(r4) - 2(F,(1-g;) + FL (1-877) Li Fo? (5.39) 


gdzie funkcje FI» FTT， Fo zależą od obszaru całkowania 


i podane są w tabelach I i II. 


Tabela 5,1 (Fir) 


min(s,t) max(s,t) 


17-Fıy 


ee nd 


۱ 2.1/2 
Powyżej ei Eet Le ) t=(1-(2,-h) ) / 


dé dla z, ۷ 2 
min(s,t) - 
ls 24 ۷ > 


(5.40) 


| Czynnik dwa w równaniu (5.59) wynika z uwzględnienia 
symetrii całkowania po przedziałach (0,d/2) i (d/2,d) 


(symetria względem punktu x4-d/? yə 


HA x dur Ay 


Pp NI - — w 
٤ , 
, ۷ 
mad > c orm > eo VEN 


< I ə. KA 
>> ار‎ ٣ - 
"i IS 7 بیو من‎ 


t 


bi ie ü Ka = عد‎ "E. 


کو 44 y‏ << ید 
" 


ELO a.‏ .ب 


Sytuacja kiedy x, e(d,2d) jest réwnowagna x, €(O,d) و‎ 


$ 
ظ‎ wystarczy tylko dokonać permutacji obszarów (permutacja ` 
x wielkości g; í Sri ) aby uzyskać szukany wynik. 


Funkcje F., F F 


Fry Fo maja swoja interpretacje fizyczna; .‏ و۲ 
przedstawiają one odpowiednio objętości części obszaru:‏ 
fazy adsorpcyjnej I , fazy adsorpcyjnej II, adsorbentu,‏ 

| które znajdują się w zasięgu oddziaływań molekuły er 

| Ilustruje to rysunek 5.7. na którym objętości te 

| zakreskowano w rćżny sposób dla dwóch możliwych | 

| przypadków położenia molekuły 1 ,.wzgledem powierzchni 


adsorbentu (płaszczyzna xy). 


| | za ELO, x 24 EL, tth) 
s WU m 
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Ryg. 5.7. Funkeje Fy, Fry: Fo jako objętości części: 
fazy adsorpcyjnej 1, fazy adsorpcyjnej II 
oraz adsorbentu,boduceğ w zasięgu oddziaływań 


molekuły 1 


Tabela 5.III podaje analityczne wzory na funkcje ' 
wystepujace w tebelach 5.1 oraz D.II z zaznaczeniem 
| przekroju objętości bryły przez te funkcje opisywanych us . 


(przekrój płaszczyzną xz). 


“abela 5.111 
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Tabela 5.111 ` szw) 


NN BENE 


L. w G(h)-G(O) 


1, = GCh)-G(2,-A1-x2) 


IA = G(zys 1-6 )-6(0) 


I, = Tx$-3x442) 


TT, 3 2 تب‎ 
Ig = 021-2221 INI 


Funkc je I,» 15, l5 podane s& w postaci róznicy 
odpowiednich wartosci funkcji G(z), będącej całką 


niecznaczong Mi cx (2- ) 2 
0 ۱ yat ls 02-4217 


G(z) = IE ax P — وس‎ 5 
A5 (2-2, a Ah- (x-X4) “_(z- 22)” 


ء x„(3-x4)‏ - —— — سب ا 


WS 


1-x4-(z-z.)” " aresin((1-x7: 


Z y - 

arcsin( Hələ) - 2x,(2-2, 
^fi-x4 

X2-2424 


"mü 
+2-2)/ 0A-x7(2-24+10)) 1 MES zl (5.41) 


Korzystağuc z tabeli 5,11 można ا یی‎ (5.39) zapisać 
w następującej postaci: 


(FTT(gT-BTT)+IC1-ET)+IO) dla 24 € (0,1) 


F(r4) ef 
2 41g t7651) *1' CI- 81) Zu € (1,n«1) 


(5.42) 
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Zgodnie z ustalonymi granicami całkowania i wzorem 


(5.33) otrzymujemy 


d/2 h h+1. 
cİN- = D PATE | F(r4)dz4)dx4 (5.43) 
O 0 h 


Uwzględniając równanie (5.42) dostaje się stąd 


d/2 h h+1 (5.44). 
۱ a 0 
2 Wy z 2(61;-67) (= KH Iden + e!Wrd 
0 O h 
gdzie h 1 h+1 


h 
2107 = (1-8; (E; Idz, + | Idz, + I'dz4) + chen + 
0 h 1 0 
1: 
+ Mgdz4 APA CZ) 
h 
Rozważając przypadek gdy x, € (d/2,d), otrzymujemy przez 


permutację Er 2 öll wyrażenie na całkę ۶۱۳ i stad 


E (5.46)‏ ليه 
I ye F.,dxjdz, + 214(478V.2)‏ 
21(W Rol = 2267-1 Beet ad AL. KB‏ 
O | `‏ 0 


Obliczymy obecnie. piervszy człon powyższego równania. 
Zgodnie z tabelą 5.1 można zapisać, że 
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Korzystając z tabeli 5.111 łatwo dowieść równości 
funkcji podcałkowych w pierwszej i drugiej całce 
występującej po prawej stronie znaku równości 
w równaniu (5.47). Zatem, 
a/ 2 t S یڈ‎ 

Pyry IX 2 G(h)dx, = G(0)dx, + aid ax, + 


0 0 0 8 
1 
ES + 1-x, )dx, l (5.48) 
t 


Całka nieoznaczona 


| | X4 
Iden E (20002... 
1- دید‎ 27 


os E 1- 2-22 "a. 
“yn aM arcsin 
MI-x7(2-24+1) afi- CEN 


" x4(arcsin( 


cəld‏ رے 
x2(x2-6) 2 df ER ey E‏ 
MAKE aros in) + NN A a +‏ 
j 1-X4 2-24‏ 
Z-Z "u Sa) =;‏ 
Azan santana (5.49)‏ )لیس 
7 


Obliczajyc przy pomocy tego wzoru całki w równaniu 
(5.48) otrzymujemy 


d/2 
4 R | h- a4 
Fr dx, = „(arcsin(h-2,) $ arcsin(z4) - — — : 
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(ena a Yo ant ab ) (5.50) 


std znajdujemy 


h d/= | (551) 
| Iren. = >(h: arcs in(h) ve A (8:9h” Lo )-75) 
00 


Nykonujac carkowanie w rövnaniu (5.45) dochodzi sie 


co następującego wzoru: 


ə a „a 3 
d = NCB) - Ster (۴ط-‎ + (61-10(77 + a ja h^ * 


le 1 
ne Pe (icm, aa), + x : (5.52) 
Analogicznie 
Wo. = Wer) (5.55) 
rt STI a 
Pocstewia jac zrównania (35.91-90.53 06 fomnia (3.5) 


a następnie 0210160 to ostatnie przez ed, otrzymuje się 


drugi współczynnik wirialny 


۱ 2. NND: ac m 

Mask ) 1+۵ - UB: + Və (2.54) 
(81-81) | ə 2 2 M. | 

(ü49h^-2h )-4-) :‏ که سس سو دس LTE (yi‏ وي UT‏ 

15 = (Wy + WIT)V2 (5.56) 


Możemy zauważyć, że 
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MEE کے‎ | (5459) 
Druzi współczynnik wirialny dla rozważanej powierzchni 
heterogenicznej (rów. 5.54) składa sie-z dwóch części: 
x niezależnej od parametru d i części us, która 
zalezy od tego parametru. Przy podstawieniu By = ETL” 
człon WA redukuje się do znanego wyrażenia na drugi 
współczynnik wirialny dla powierzchni homogenicznej 
/*^,292/. Człon W równa Się wówczas zeru. W opisuje 
zatem oddziaływania molekuł znajdujacych się na różnych 
płatach powierzchni, Gdy d— co Wówczas wartość tego 
członu maleje do zera, Odpowiada to fizycznie 
powiększaniu płatów powierzchni do momentu kiedy 
oddziaływania między molekułami zaadsorbowanymi na 
różnych płatach stajy się zaniedbywalne. Otrzymaliśmy 

w ten sposób sformułowanie modelu piatowezo powierz- 
ohni heteroxeniczned. Wyniku stud,ze grupowanie w piaty 
centrów adsorpcy jnych jest mozliwe jeśli 32 
między centrami scharakteryZowanyni róznymi enerziami 
adsorpcji su duże. Dla rzeczywistych adsorbentow 
odległości te są równe w przybliżeniu zakresowi. 
oddziaływać międzycząsteczkowych, posłużziwanie sie więc 
(szczególnie przy badaniu adsorpcji z bocznymi oddziały- 
waniami) modelem płatowym powierzchni wymuga ostrożności. 
Z równania (5.55) wynika, że wartość współczynnika Ha 
zależy od مرن‎ , można więc przypuszczać, zek 
zarówno wartości potencyaiow Vr, ۷ jak 1 tenpera- 


tura będę miały na nig duzy wpływ: 
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Tabela. IV 


— — = “Un == — 5 
L ———— — —— — s ———— —— ہے‎ ee | 


| 
układ | T = 85 K 100 K 120 K 150 K 200 K 250 K 


= İ 
| 一 一 一 一 一 — 


m" 263448 | 495446 


WS 1296447 | 466045 
Vi 52250 
% |12.68 | 10.38 
AS š 
-W5 | 486446 
864548 


21.6 


877449 | 182+13 
0ءع ۷ 

359+12| 1716+9 
Vi, 220 

24.55. 24,3 
سج‎ m 


۷ tabeli IV  podane sy wartoáci współczynników Wo 1 Wo 
obliczone na podstawie równań (5.52 - 5.56) oraz 
procentowy wkiad və do W. ) وا‎ 2832). Obliezenia 
przeprowadzono dla trzech układów określonych parami 
potencjałów (vj, Vij) i sześciu wartości temperatury 
zamieszczonych w tabeli, przyjmując następujące wartości 
parametrów: 422 i h-1. 

Zgodnie z oczekiwaniami, procentowy udział ws W Wo 


rośnie ze wzrostem różnicy potencjałów Vn Vinge aer male je 


ze wzrostem tem;eravury.Daje sie przy tym zauważyć, że 


wpływ temperatur; na ten udzial jest gest niewielki 
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przy dostatecznie dużej różnicy potencjałów. Widoczna 
jest również w tym przypadku zbieżność wkładu procento- 


wego wynikajuca z następującego oszacowania: 


yt 288117 STARA „028 


Wo 


= (4.047 + y 4 0.24? (5.60) 


(&i- RL, 


Nierówność (5.60) została otrzymana z równań (5.54, 
5,55) dla wartości 0-2 i hal 

W5 nigdy nie przekracza więc wartości 0.247 4H] ) lub 

0.359 ( . Y sytuacji kiedy ولا‎ jest zaniedbywalne, 

“2 jest tym bardziej bez znaczenia. 


Rozważmy obecnie izotermę adsorpcji: 
BE Kun + 2W_n* (5.61) 
H 2 SÉ 


Stara Henryego dla rozpatrywanej dwucentrowej powierz- 


chni heterozenicznej wynosi: 


hd 20 
1 
Ku = > (e(r3-1)dr )مس اج‎ 3 Gu 1(0 + | rr Daaa 5 
Kä 00 d 
= 0.5*h(gp *£11-2) (5.62) 
Można tez ja zapisać w postaci sumy: 


gdzie K+ i Ky; SA stazymi Henry'ego dla powierzchni 
homogenicznych z potenc jatami odpowiednio vr 1 ۰۶ 
Korzystajuc z wyprowadzonego wyrażenia (5.63) i równania 


(5.54) na drugi współczynnik wirialny ۷2 otrzymujemy 
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następującą postać izotermy adsorpcji 


RZĘS: 


d 02 o „2 
> Kin + en 5 br: = Kyyn + 2 


(5.64) 


(5.65) 


(5.66) 


(5.67) 


Powyżej, wielkości RE i izə sy jak łatwo można 


zauwazyó, izotermami adsorpcji dla powierzchni homo- 


genicznych . Zatem człon ۳ 


opisuje adsorpcję na 


powierzchni heterogenicznej przy założeniu płatowej 


topografii powierzchni,natomiast [' jest dodatkową 


wartością adsorpcji wynikejuca z uwzględnienia oddzia- 


ływań między molekułami zaadsorbowanymi na miejscach 


adsorpcyjnych o różnej energii adsorpcji. Oczywiście 


wartość J’ związana jest ze współczynnikiem m 


Fakt, , ze (H wyraża się jako średnia arytmetyczna. 


stężeń powierzchniowych 2 1 Kë jest niczym innym, 


jak odzwierciedleniem równości szerokości pasów - d - 


obszarów scharakteryzowanych przez potencjał oddziaiy- 


wania gaz-ciało stałe v(r) 


. € przypadku 


¿ay na 


powierzchni adsorbentu mozna wyróżnić wiele takich 


obszarów , izoterma adsorpcji wyraża się analogicznym 


do (5.64) równaniem: 


| `- Że; [3 z Ké Dä 


(5.65) 
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gdzie sumowanie przebiega po wszystkich vyröznionych 
przez nas obszarach "i". Powyżej m jest stezeniem 
powierzchniowym dla homogenicznego płatu, określonego 
przez potencjał NĄ: 207 dest wartością adsorpcji. 
wynikajaca z uwzględnienia kontaktu. między. molekułami 
zaadsorbowanymi na i-tym i  j-tym płacie powierzchni, 
c, jest związane z rozkładem potencjału v(r) na 
powierzchni. Fierwszy człon równania (5.65) jest 
powszechnie znanym wyrażeniem opisującym adsorpcję 

na płatowo-heterogenicznej powierzchni, ściśle określo- 
nym przez parametry potencjału wyróżnionych płatów, 
natomiast obliczenie członu drugiego wymaga dodatkowo 
informacji o ich położeniu i kształcie. Znajomość 
topografii powierzchni umożliwia więc obliczenie 
wkładu Î ; którego wielkość jest jednocześnie miarą 


odchyleń topografii od skrajnego, pzatowego modelu 


powierzchni niejednorodnej. Zatem rozmieszczenie centrów ` 


adsorpcyjnych, przy którym ۲ osiega wartość maksy- 
malny może być uważane za przypadek topografii "random". 
Badania M. Jaroíca i M. Borówko 24/ wykazały, ze izo- 
termy adsorpcji dla dwucentrowej powierzchni heteroge- 
nicznej zawarte sy między dwiema skrajnymi izotermami 
odpowiadającymi pewnym "ekstremalnym" rozmieszczeniom 
centröw PER" jednym Es Mer R MUR 1 
płatowy. 

Ilustracyjne obliczenia numeryczne do tego rozdziału 


przeprowadzono dla układów scharakteryzowanych przez 


potencjały: v4-2000 cal/mol, 0ء ۔۔۷‎ cal/mol  oraz 


DU g " 
prt > 


-— = =. 
N 1 LAS u d 
k. Y jus کے‎ i 
| B T | 


3 ا 


۷۰۶۵000 cal/mol, Vip” 2500 cal/mol, używając wartości 
współczynników wirielnych zamieszczonych w tabeli IV 
dla^dwóch wartości temperatur T-125 K i Ta 50 K. 
Rysunek 5.8.4 przedstawia izotermy adsorpcji obli- 
czone zgodnie z równaniem (5.61) w porównaniu 

z 120 terms Henry'ego (Xun) dla powierzchni heteroze- 
nicznej, rys.5.8.b przedstawia wartości | + (rów.5.67), 
rys. 5.8.c pokazuje procentowy udział | ۱ .100% i 
Zgodnie z dysrusig dotycząca 0 współczynnika 
wirialnego (str. 125, 126),["* przyjmuje duże WAREN 
w przypadku gdy różnica potencjarów Vy i Vir jest 
duża. Zależność | + NB od temperatury i gęstości 
adsorbatu (ciśnienia) jest analogiczna. 

Podkreálmy tu raz jeszcze, że zaniedbanie pe di rea 
(5.64) jest równoważne z przyjęciem modelu płatowego 
tovozrafij powierzchni (molekuły zaadsorbowane na 
centrach o różnej energii adsorpcji aie oddziatrywuj& 

ze S0bu). Wartość |" ‘Cw dyskutowanym przypadku - 


adsorpcji sztywnych kul) jest zawsze dodatnia, zatem 


adsorpcja na powierzchni wykazujęcy pośrednią topografię 


miejsc adsorpcyjnych Jest większa niż adsorpcja na 
powierzchni. z topografin piatową. Procentowy udział 
| ۱ ۵ maleje ze wzrostem temperatury zgodnie 
z Tövnaniem : 


d i 3 
Limes (exp(vy/xT) - exp(v,/ kT)) = 0 (5.66) 


Przy dość małej różnicy potencjałów Verve GPO SEES 


ehnia prawie homogenicana) , odpowiadającej szerokości 
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rozkładu energii adsorpcji dla układu argon-ruty1/ 9% , 
udział | * molekuł zaadsorbowanych na centrach o różnej 
energii adsorpcji wynosi kilka procent. Należy się więc 
spodziewać, że w przypadku pad SRA rzeczywistych o dużej 
niejednorodności i dużej średniej energii ,adsorpcja może 
być znacznie większa. W tych szczególnych warunkach | * 


(topografia powierzchni) nie może być zaniedbywana. 
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PODSUMOWANIE 
Rozdział I: 


W rozdziale I przedstawiono dotychczasowy stan badań 
nad efektami topografii powierzchni heteregenicznych. 
zaskakujący jest fakt, że prawie 90% prac dotyczących 
teorii adsorpcji na powierzchniach neterogenicznych 
bazuje na płatowym modelu powierzchni, podczas gdy 
ścisłe określenie własności układów z heterogeniczną 
powierzchnię wymaga kompletnej znajomości topografii 
całej powierzchni. W tej sytuacji stało się konieczne 


systematyczne przebadanie tego zagadnienia, 
Rozdział II: 


Do opisu topografii powierzchni wprowadzono funkcję 
lokalnego, różniczkowego rozkładu energii adsorpcji 
X (8, رظ‎ - różniczkowy rozkład energii adsorpcji E, 
na heterogenicznym płacie najbliższym sąsiedztwie 
centrum adsorpcyjnego o energii E. 
Dla idealnie płatowej topografii powierzchni hetero- 
geniczne ل‎ 

X, (EE) = (E, EI 
natomiast w przypadku całkowicie przypadkowego 
rozkładu miejsc adsorpcy jnych, X, (2, رظ‎ = XE.) 
Porównano wielkoáci termodynamiczne obliczone dla 


trzech typów topografii powierzchni heterogenicznych 
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ptatowej, pośredniej i "random" przy użyciu przybliżenia 


Bragga-Williamsa. 


Rozdział III 


metoda pezvolitania. kondensacy metro 95 ۳ ۵‏ دكن اتن رأ ند 
tyczniə słusznego przybliżenia kondeusaey inego‏ 

۷۸۷۷ 77:207 0120 różniczkowe rozkłady energii adsorpcji 

dla powierzchni heterozeniczneğ typu "random" i płatowej. 


kozdziaz IV 
Uzyvağuc przyblizenia quasi-chemicznezo FG, przeprowa- 
dzono systematyczne badania adsorpcji dla powierzchni 
heterosenicanej "random" w obszarze temperatur: 


superkrytycznym i subkrytycznyn. 
Rozdział V 


Sformułowano trójwymiarowy model adsorpcji mobilnej 

na powierzchni nie jednorcanej. Wyprowadzono równania 
na wyznaczanie gęstości lokalnej Da dla powierz- 
chni heterogeniczneğ oraz przedyskutowano zbieżność 
całkowych vövnaü adsorpcji dla modelu dwu i trój- 
wymiarowego. wyprowadzono drugi współczynnik wirialny 


dla dwucentrowe j powierzchni heterogenicznej. 
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W przyblizeniu Brasga-Williamsa, informacja 
i o topografii powierzchni potrzebna jest do obliczenia 
potencjału średniej siły działającej na zaadsorbowaną 
molekułę. Wystarcza do tego znajomość funkcji rozkładu : 
Xj(E ,E) przy pomocy której oblicza sie wzgledne 
pokrycie obszaru-najbliższego sąsiedztwa miejsca 
adsorpcyjnego o energii E . Jak się okazuje, przyjęcie 
na X, (E, E) globalnej funkcji rozktadu energii 
adsorpcji X(E) prowadzi do wyższych wartości izoterm 
adsorpcji i ciepeł izosterycznych niż założenie, że 
£= ó(E,-E). Parametrem "okreélajacym" topografie 
powierzchni jest w tym przypadku szerokość D funkcji 
rozkładu X,(E,,E). Gdy D = 0 wówczas X, = Ś(E,-B) 
i mamy do czynienia z przypadkiem powierzchni płatowė j. 
Małe wartości D oznaczają istnienie obok siebie 
( w płatach) miejsc adsorpcyjnych o zbliżonej energii 
adsorpcji. W miarę wzrostu wartości D od O’ do 00 
(do szerokości funkcji rozkładu energii X(E)) przecha- . ` 
dzimy od topografii płatowej przez pośrednią do random. 
Wpływ parametru D na wartość adsorpcji i izosteryczne 
ciepła adsorpcji przejawia się w ich wzroście przy 
| zwiększeniu wartości D. 
| Posługując się modelem trójwymiarowym adsorpcji mobilnej, 
| możemy modelować powierzchnię heterogeniczną wprowadzając 
z góry zadaną funkcję potencjału v(x,y,2). Zależy ona "s 
najczęściej od pewnych parametrów określających topo- 
grafię powierzchni (kształt potencjału VUX NAE 


| Parametry te, decydują o przebiegu potencjału Med 
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a w konsekwencji o przebiegu gesto6oi l oxelne. KARES 
od której zależ4 vielkoöci termodynamiczne Taxim 


perametrem topografii jest w.naszym przypadzu "a", 
Parametr ten ma ściśle określony sens fizyczny staiej- ° 
sieciowej. Zwiekszeniu wartoSci "a" towarzyszy zmnie j- 
szenie nie jednorodności lokalnej -obszarów o okreálonym 
promieniu (zasięgu). Sytuacja taka odpowiada więc, 
tworzeniu się płatów złożonych z miejsc 20150252037 

o zbli&one; eESTTI ad 02p2J1 (relucsa DRO ET 
dziale II ). W miarę zmniejszania sie wartości para- 
metru "a" mamy do czynienia ze zwiększaniem niejedno- 
rodności, a więc proces odwrotny do poprzedniego (ana- 
logiczn'" ¿o  DL—»00). Możemy w tym przypadku mówić 

o przechodzeniu do toşografii powierzchni typu "random". 
Badağuc efekt topografii na adsorpcję, obserwujemy 
analogiczny wpływ jak poprzednio tj. ze wzrostem wartości 
parametru topografii "a" otrzymuje sie zmniejszenie 
ilości zaadsorbowanej. Zauwaza sie również, ze wpływ. ten 
jest większy w przypadku dużych wartości średniej 

energii 2 ° 

Model trójwymiarowy dostarcza także tłumaczenia tego 
faktu posługując się gęstościy lokalną. mianowicie, 
można stwierdzić, ze ze wzrostem wartości varametru "a" 
maleje |liczba czystek zaadsorbowanych na miejscach 
najbliżej suSiadujycym z miejscem wyróżnionym. Tym e 
samym maleje udział adsorpcji wynikający 2 odds Tatty wit 
miedzycze,steczkowych w fazie adsorpcyjnej. 

W celu ilościowegc zbadanie zależności (wzrostu czy 
obniżenia adsorpcji) od średniej energii adsorpcji, 
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niejednorodności i stałej sieciowej przebadano 
adsorpcję na dwucentrowe j powierzchni heterogeniczne j 

o periodycznym potencjale określonym przez stałą 
sieciową d. Obliczenia wykazały zgodność z poprzed- 
nimi wynikami. Stwierdzono jednocześnie, ze odchylenie 
topografii od modelu piatowego scharakteryzowane jest 
przez wielkość oddziatywaui cząsteczek zaadsorbowanych 
na miejscach adsorpcy jnych o różnej energii adsorpcji, 
która zależy bezpośrednio od .otencjału gaz-ciało stałe 


i parametru topografii d. 


Badania przeprowadzone w trójwymiarowym modelu 
adsorpcji dotyczy najprostszego przypadku adsorpcji 
sztywnych kul. [Interesujące jest przeprowadzenie 
analogicznych badań dla bardziej sxomplikowanych 
(uwzsględniajycych przyciąganie między cząsteczkami) 
potencjałów  gaz-gaz. Badania takie pozwoliiyby 
otrzymać informacje a subtelnych efektach takich 


jak np. przejścia fazowe. 
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SPIS WAZNIEJSZYCH SYMBOLI 


À - energia swobodna Helmholtza 

a - stała sieciowa (parametr topografii) 

Cy - pojemność cieplna 

6 - liczba najbliżej sąsiadujących miejsc adsorpcy jnych 


c - indeks, oznaczający wartość w punkcie krytycznym 


D - szerokość funkcji X4 (parametr topografii) 

d - szerokość pasów (parametr topografii) 

E - enerzia adsorpcji 

E - energia wewnętrzna rozdz. ۱ 3. 

p - i-czysteczkowa funkcja rozdziału Ursella 

Le - funkcja rozdziału dla wewnętrznych stopni swobody 
molekuły gazowej 

f, - funkcja rozdziału dla zaadsorbowane j molekuły 


f(x, ,) - funkcja Mayera 

۽٣2) (ر‎ - radialna funkcja dystrybucji 

g(r, no) - dwucząsteczkowa funkcja korelacji 

h - stała Planca 

h - parametr określający zakres potencjału v(r) 
rozdz. 5.5 (parsmetr topografii) | 
K - stała Langmuira 


K.. - stała Henry ego 


H 
k - stala Boltzmana 

M - ilość miejsc adsorpcy jnych 

m - masa cząsteczki 

N - ilość zaadsorbowana na powierzchni 


n - gęstość adsorbatu 
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(ri), n CE; 3? ~ jedno i dwucząsteczkowa 
funkcja rozkładu 


p - ciśnienie adsorbatu 


R - stała gazowa 


r - parametr niejednorodności 
r - wektor położenia 
- entropia rozdz. IV.5 


S 
S - pole powierzchni adsorbentu 
T - R مان‎ thie 

T 


e 

U - extra energia(Cerofoliniego) 

NT energia potencjalna N cząsteczek 

2- drugi współczynnik wirialny dla powierzchni, 

heterogeniczneğ 

W5 - drugi współczynnik wirialny dla powierzchni 
heterogenicznej o topografii płatowej 

Vİ - człon do W. wynikający z oddziaływń dwóch ہے‎ 
molekuł zaadsorbowanych na centrach o różnej 
energii adsorpcji 

V - względne pokrycie powie zar 

V - objętość układu adsorpcyjnego rozdz. V 


V 


os? 


nes granice obszaru kondensacy jnego 


o v(r) - potencjał oddziaływania gaz-ciało stałe 


W - energia oddziaływania dwóch molekuł zaadsorbowanych 
. na miejsaach najblizej susiadudacych ze Sobą 

u) - potencjał oddziaływania gaz - gaz 

X(E) - funkcja rozkładu energii adsorpcji 

TE - ciepto izosteryczne 


Zy - kanoniczna suma stanów 


Toe» Tze ~ dwu i trzywymiarowa temperatura krytyczna 
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€ - parametr określający potencjał v (z) - średnia 
7 energia adsorpcji dla powierzchni z potencjałem 71 
Ç - parametr określający wielkość potencjatu 
| 2 20 ۱ 0 
(Jp, E) - równanie izotermy: adsorpcji na miejscach 
adsorpcy jnych scharakteryzowanych energią E 
Par potencjał chemiczny składnika w fazie adsorpcy jnej 
ə potencjał chemiczny gazowego adsorbatu ۰ 
C - średnica cząstki adsorbatu 
| m - nadmiar powierzchniowy Gibbsa 
| | + - wkiad do ۴ wynikający z uwzzlędnienia oddziaływań 
cząsteczek zaadsorbowanych na różnych miejscach 
adsorpcyjnych (o Mone enerzii adsorpcji) 
^C - izotermiczny współczynnik ściśliwości 
d - funkcja zeta Riemanna 
N- przedział zmienności wartości energii adsorpcji E 


Tt- stała matematyczna 
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